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Abstract 

We consider the motion of a particle governed by a weakly random Hamiltonian flow. We 
identify temporal and spatial scales on which the particle trajectory converges to a spatial Brown- 
ed ' ian motion. The main technical issue in the proof is to obtain error estimates for the convergence 
of the solution of the stochastic acceleration problem to a momentum diffusion. We also apply 
our results to the system of random geometric acoustics equations and show that the energy 
, density of the acoustic waves undergoes a spatial diffusion. 
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The long time, large distance behavior of a massive particle in a weakly random time-independent 
potential field is described by the momentum diffusion: the particle momentum undergoes the 
Brownian motion on the energy sphere. This intuitive result has been first proved in j§] in dimensions 
higher than two, and later extended to two dimensions with the Poisson distribution of scatterers in 
[2j . On the other hand, the long time limit of a momentum diffusion is the standard spatial Brownian 
motion. Hence, a natural question arises if it is possible to obtain such a Brownian motion directly 
as the limiting description in the original problem of a particle in a quenched random potential. This 
necessitates the control of the particle behavior over times longer than those when the momentum 
diffusion holds. 

We consider a particle that moves in an isotropic weakly random Hamiltonian flow with the 



Or 



X 



Hamiltonian of the form H$(x, k) = Ho(k) + ySHi(x, k), k = |k|, and x, k £ WL d with d > 3: 
dX^ rlK^ 

- 1F = y^H s , — = -V x H s , X 5 (0) = 0, K 5 (0)=k . (1.1) 

Here Ho(k) is the background Hamiltonian and iii(x, fc) is a random perturbation. As we have 

k 2 r 

mentioned, it has been shown in [S] that, when ^(x, k) = — + v SV (x), and under certain mixing 

assumptions on the random potential V(x), the momentum process K s (t/8) converges to a diffusion 
process K(t) on the sphere k = ko and the rescaled spatial component X s (t) = 5X s (t/5) converges 
to X(t) = Jq K(s)ds. This is the momentum diffusion mentioned above. Another special case, 

ff{(x,k) = (cD+VMx))|k|, (1.2) 
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arises in the geometrical optics limit of wave propagation. Here Co is the background sound speed, 
and ci(x) is a random perturbation. This case has been considered in pQ, where it has been shown 
that, once again, K s (t/5) converges to a diffusion process K(t) on the sphere {k = ko} while 
X s (t) = SX 5 {t/5) converges to X(t) = Co J* K(s)ds, K(t) := K{t)/\K{t)\. 

We show in this paper that this analysis may be pushed beyond the time of the momentum 
diffusion, and that under certain assumptions concerning the mixing properties of H\ in the spatial 
variable there exists ao > so that the process converges to the standard Brownian 

motion in R d for all a € (0, ao). The main difficulty of the proof is to obtain error estimates in the 
convergence of K s (-) to the momentum diffusion on time scales of the order 5~ 1 . The error estimates 
allow us to push the analysis to times much longer than <5 _1 where the momentum diffusion converges 
to the standard Brownian motion. The method of the proof is a modification of the cut-off technique 
used in and jH]. 

A similar question arises in the semi-classical limit of the quantum mechanics and high frequency 
wave propagation. The Wigner transform [""], or the phase space energy density of the solution of the 
Schrodinger equation, is approximated in a weakly random medium by the solution of a deterministic 
linear Boltzmann equation [3]. This behavior is also conjectured for the acoustic and other waves 
in a weakly random medium [""""". As in the momentum diffusion model for a particle, the long time 
limit of the Boltzmann equation is the spatial diffusion equation. It has been recently shown in j3J 
that, indeed, one may push the analysis of [3| beyond the times on which the Boltzmann equation 
holds and obtain the diffusive behavior of the energy density of the solutions of the Schrodinger 
equation in the weak coupling limit. 

We also apply our results to the problem of multiple scattering of the acoustic waves. Our 
approach is different from that of [3] mentioned above: we first consider the random geometrical 
optics approximation of the wave phase space energy density. The rays in the phase space satisfy 
the Hamiltonian equations "l.ljl with the Hamiltonian given by (|1.2|) . Therefore, the aforementioned 
convergence result of the solutions of (jl. 1|) to the standard Brownian motion, combined with the 
error estimates on the geometrical optics approximation of the Wigner distribution of the solutions 
of the wave equation, allows us to establish rigorously the diffusive behavior of the wave energy 
density. To the best of our knowledge, this is the first result of such kind for classical waves. 

This paper is organized as follows. Section [21 contains the main results on the convergence 
of solutions of (jl.lj) to the Brownian motion as well as the assumptions on the random medium. 
Sections 01 and 0] contain the proof of our main results, Theorems 12 . 1| 12*31 and [2*771 the error estimates 
on the passage to the momentum diffusion and the passage from the momentum diffusion to the 
spatial diffusion, respectively. Section El discusses the application to the wave equation. Finally, 
Appendix [J3 contains the proof of a technical result that is stated in Lemma 13.51 We note that all 
constants appearing throughout the paper do not depend on 5 £ (0, 1] unless otherwise specified. 

Acknowledgment. The research of TK was partially supported by KBN grant 2PO3A03123. 
The work of LR was partially supported by an NSF grant DMS-0203537, an ONR grant N00014-02- 
1-0089 and an Alfred P. Sloan Fellowship. 

2 The main result and preliminaries 
2.1 The notation 

As we will avoid the sine ular point k = 0, we denote R d := R d \ {0} and M% d ■= R d x R d . Also 
§^ _1 (x) (Br(x)) shall stand for a sphere (open ball) in R d of radius R > centered at x. We shall 
drop writing either x, or R in the notation of the sphere (ball) in the particular cases when either 
x = 0, or R = 1. For a fixed M > we define the spherical shell A(M) := [keif: M" 1 < |k| < M] 
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in the k-space, and A(M) := M. d x A(M) in the whole phase space. Given a vector v G Rf we denote 
by v := v/|v| G S^" 1 the unit vector in the direction of v. For any set A we shall denote by A c its 
complement. 

For any non-negative integers p,q,r, positive times T > T* > and a function G : [T*, T] x — > 
R that has p, q and r derivatives in the respective variables we define 

W G \&--=J2 SU P |^fl£G(*,x,k)|. (2.1) 

(t,x,k)e[T„T]xM M 

The summation range covers all integers < a < p and all integer valued multi- indices \/3\ < q 
and |7| < r. In the special case when T* = 0, T = +00 we write ||G|| P)9ir = HGHp^- 00 ^. We denote 
by Cf ,q,r ([0, +00) x M% d ) the space of all functions G with ||G|| P)gjr < +00. We shall also consider 
spaces of bounded and a suitable number of times continuously differentiable functions C^ ,q ( 
and C%(M%) with the respective norms || • \\ Pt q and || • || p . 



2.2 The background Hamiltonian 

We assume that the background Hamiltonian Ho{k) is isotropic, that is, it depends only on k = |k|, 
and is uniform in space. Moreover, we assume that Hq : [0, +00) — > K is a strictly increasing function 
satisfying Hq(Q) > and such that it is of C 3 -class of regularity in (0, +00) with H' (k) > for all 
k > 0, and let 

h*(M):= max (H' (k) + \H'^k)\ + \H'^'(k)\), K(M) := mm H' (k). (2.2) 

Two examples of such Hamiltonians are the quantum Hamiltonian H$(k) = k 2 /2 and the acoustic 
wave Hamiltonian Ho(k) = c$k. 



2.3 The random medium 

Let (fi, £, P) be a probability space, and let E denote the expectation with respect to P. We 
denote by ||X||£j>(q) the L p -norm of a given random variable X : Q, — > M, p G [1, +00]. Let iJi : 
M rf x [0, +cxo) x Q —>■ M be a random field that is measurable and strictly stationary in the first variable. 
This means that for any shift x € k G [0, +00), and a collection of points xi, . . . , x n G R rf the laws 
of (iJi(xi + x, fc), . . . , i?i(x n + x, /c)) and (-ffi(xi, fc), . . . , i?i(x n , fc)) are identical. In addition, we 
assume that Ei?i(x, fc) = for all A; > 0, x G W 1 , the realizations of i?i(x, A;) are P-a.s. C 2 -smooth 
in (x, k) £R d x (0, +00) and they satisfy 

Dij(M) := max ess-sup |fl£8j[.Hi(x, < +00, i,j = 0,1,2. (2.3) 

l Q l=* (x,fc,o;)eM [i x[Af- 1 ,Af]xn 

We define D(M) := Eo< l+i <2 Aj(M). 

We suppose further that the random field is strongly mixing in the uniform sense. More precisely, 
for any R > we let and CJj be the cr-algebras generated by random variables H\(x,k) with 
k G [0, +00), x G Mr and x G M R respectively. The uniform mixing coefficient between the u-algebras 
is 

cj>(p) := sup[ |P(fl) - F(B\A)\ : R > 0, A G C R , B G C| +p ], 
for all ^ > 0. We suppose that cfr(p) decays faster than any power: for each p > 

hp := sup/AKp) < +00. (2.4) 

p>0 
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The two-point spatial correlation function of the random field Hi is R(y, k) := E[i?i(y, k)H\(0, k)}. 
Note that (|2.4|) implies that for each p > 

4 

hp(M) :=VV sup (1 + \y\ 2 f /2 \d$R(y, k)\ < +oo, M > 0. (2.5) 

i=0 | a |=i (y,fc)6K d x[A/-i,A/] 

We also assume that the correlation function i?(y, Z) is of the C°°-class for a fixed Z > 0, is sufficiently 
smooth in Z, and that for any fixed / > 

R(k, Z) does not vanish identically on any hyperplane H p = {k : (k • p) = 0}. (2-6) 

Here i?(k, Z) = J -R(x, Z) exp(— ik • x)cZx is the power spectrum of 

The above assumptions are satisfied, for example, if H±(x,k) = ci(x)h(k), where ci(x) is a 
stationary uniformly mixing random field with a smooth correlation function, and h(k) is a smooth 
deterministic function. 



2.4 Certain path-spaces 



For fixed integers rf,m > 1 we let C 



d,m ._ 



C([0, +oo); M. d x K™): we shall omit the subscripts in the 



as the canonical 



notation of the path space if m = d. We define (X(t),K(t)) : C 
mapping (X(t; ir), K(t; %)) := vr(t), ir G C d,m and also let s (tt)(-) := %(■ + s) be the standard shift 
transformation . 

For any u < v denote by M v u the cr-algebra of subsets of C generated by (X(t),K(t)), t £ [u,v]. 
We write Ai v := A^q and A4 for the a algebra of Borel subsets of C. It coincides with the smallest 
cr-algebra that contains all M*, t > 0. 

Let 6*(M) := H Q (M' 1 ) /(2D(M)). For a given M > and 5 G (0, <5*(M)] we let 



M 6 := max <| H Q l {H {M) + 2v / 5L»(M)), i^o" 1 ( ff ( ^ ) - 2VdD(M) 



(2.7) 



For a particle that is governed by the Hamiltonian flow generated by H$(x., k) we have M 5 < 
1-^(^)1 < Mg for all t provided that K(0) G yl(iW). Accordingly, we define C(T, 5) as the set of paths 
vr G C so that both (2M5)" 1 < |Jf(i)| < 2M 5 , and 



X(t) - X(u) - / H' (K(s))k(s)ds 



< D(2M s )\/5(t - u), for all < u < t < T. 



In the case when 5 = 1, or T = +oo we shall write simply C(T), or C(5) respectively. 



2.5 The main results 

Let the function (j>g(t,x., k) satisfy the Liouville equation 

^ + V X F 5 (x, k) • Vk^ - V k H s (x, k) • V x / = 0, (2.8) 
/(0,x,k) =<M<*x,k). 

We assume that the initial data ^o( x > k) is a compactly supported function four times differentiable 
in k, twice differentiable in x whose support is contained inside a spherical shell A(M) = {(x, k) : 
M' 1 < |k| < M} for some positive M > 0. 
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Let us define the diffusion matrix D mn by 

r, n 1 f°° d 2 R(H'(l)sk,l) , 1 Z" 00 d 2 R(sk,l) , , . . 

D mn (k,l) = -- A ow ^ dg = _^^ / v m,n = l,...,d. 2.9 

2 y_oo ox n dx m 2H' (l) J ^ dx n dx m 

Then we have the following result. 

Theorem 2.1 Zei a 6e the solution of and let <f> satisfy 



E ^(D mn (kk)^+H' (k)k-V x $ (2.10) 
0(0, x,k) = 0o(x,k). 



9t 

m,n=l 



Suppose that M > Mq > and T > Tq > 0. Then, there exist two constants C, ao > swc/i i/iai for 
all T >T 



sup 

(t,x,k)6[0,T]xK 
for all compact sets K C A{M). 



<CT(l + ||0 o || M )5 ao (2.11) 



Remark 2.2 We shall denote by C, Ci, . . ., ao, «i, . . ., 7o, 7i, • • • throughout this article generic 
positive constants. Unless specified otherwise the constants denoted this way shall depend neither 
on 5, nor on T . We will also assume that T > Tq > and M > Mq > 0. 

Remark 2.3 Classical results of the theory of stochastic differential equations, see e.g. Theorem 6 
of Chapter 2, p. 176 and Corollary 4 of Chapter 3, p. 303 of H, imply that there exists a unique 
solution to the Cauchy problem Ij2.10|) that belongs to the class C 6 ' ' ([0, +oo) x R^). This solution 
admits a probabilistic representation using the law of a time homogeneous diffusion Xj k whose 
Kolmogorov equation is given by (|2.10|) . see Section l3~Kl below. 

Note that 

^ D nm{ k, k)k m - ^ J dxnQxm k m ds 



m=l m 



i_ r°° d_ ( dR(sk,k) \ 

^HWLndsy d Xn ) S ~ 

and thus the -fT-process generated by (|2.10|) is indeed a diffusion process on a sphere k = const, or, 
equivalently, equations (|2.10j) for different values of k are decoupled. Assumption (|2.6|) implies the 
following. 

Proposition 2.4 The matrix D(k, I) has rank d — 1 for each k 6 § d_1 and each I > 0. 

The proof is the same as that of Proposition 4.3 in pQ. It can be shown, using the argument given 
on pp. 122-123 of ibid., that, under assumption l|2.6j) . equation 1)2.10(1 is hypoelliptic on the manifold 
R d x S^ 1 for each k> 0. 

We also show that solutions of (|2.10j) converge in the long time limit to the solutions of the spatial 
diffusion equation. More, precisely, we have the following result. Let 7 (t,x, k) = ^(t/j 2 , x/7, k), 
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where <fi satisfies (|2.1()jl with an initial data <^ 7 (0, t, x, k) = </>o(7x, k). We also let w(t,x.,k) be the 
solution of the spatial diffusion equation: 

m,n=l " 171 

io(0,x,fc) = ^ (x, k) 

with the averaged initial data 

&(x,A0 = -!- / <fo(x,k)dfi(k). 

Here cZQ(k) is the surface measure on the unit sphere S^ -1 and T n is the area of an n-dimensional 
sphere. The diffusion matrix A := [a nm ] in (|2.12|) is given explicitly as 

a nm (k) = f H' (k)k nX m(k)dn(k). (2.13) 

The functions Xj appearing above are the mean-zero solutions of 

E J- (^&*>ir) = -^wv ( 2 - 14 ) 

m,n=l ^ 



Note that equations ()2.14|) for Xm are elliptic on each sphere {|k| = k}. This follows from the fact that 
the equations for each such sphere are all decoupled and Proposition 12.41 Also note that the matrix 
A is positive definite. Indeed, let c = (c l5 . . . , c^) E M. d be a fixed vector and let Xc '■= El=i c mXm- 
Since the matrix D is non-negative we have 

(.4c,cw = -T^jtj^ *c(M)^ (o-CtO^S) «*) (2-15) 

^EX-<M^M^W-o* 



7§d-l 



m,n=l ' 

(D(k,OV X c(k,0,Vxc(k,0)R d ^(k) > 0. 



'Idr-X 



The last equality holds after integration by parts because -D(k, Z)k = 0. Moreover, the inequality 
appearing in the last line of ()2,15j) is strict. This can be seen as follows. Since the null-space of the 
matrix D(k, I) is one-dimensional and consists of the vectors parallel to k, in order for (Ac, c) K <! to 
vanish one needs that the gradient Vx c (k, I) is parallel to k for all k S § rf_1 . This, however, together 
with (|2.14j) would imply that k • c = for all k, which is impossible. 
The following theorem holds. 

Theorem 2.5 For every < T* < T < +oo the re-scaled solution 7 (t,x, k) = (f>(t/^/ 2 , x/7, k) of 
converges as 7 — > in C([T*, T]; L°°(M. 2d )) to w(t, x, k). Moreover, there exists a constant 
C > so that we have 

\\w(t, ■) - fa(t, 0110,0 < C ( 7 T + ^7) l|0o||i,i (2-16) 

/or a// <t <T . 
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Remark 2.6 In fact, as it will become apparent in the course of the proof, we have a stronger result, 
namely T* can be made to vanish as 7 — > 0. For instance, we can choose T* = 7 3 ^ 2 , see (|4.16|) . 

The proof of Theorem 12,51 is based on some classical asymptotic expansions and is quite straightfor- 
ward. As an immediate corollary of Theorems 12.11 and 12.51 we obtain the following result, which is 
the main result of this paper. 



Theorem 2.7 Let fa. be solution of (\2.8\i with the initial data fa(0, x, k) = <^o(<5 1+Qx ) k) and let 
w(t,x) be the solution of the diffusion equation $2.1<$f with the initial data w(0,x,k) = fa(x.,k). 
Then, there exists ao > and a constant C > so that for all < a < ao and all < T* < T we 
have for all compact sets K C A(M) : 



sup 



\w(t,x,k) -Efa(t,x,k)\ < CT5 a °- 



(2.17) 



where fa(t,x,k) := fa (t/5 1+2a , x/5 1+a , k) 



Theorem 12.71 shows that the movement of a particle in a weakly random quenched Hamiltonian is, 
indeed, approximated by a Brownian motion in the long time-large space limit, at least for times 
T <C 5~ a ° . In fact, according to Remark l2.6l we can allow T* to vanish as 5 — ► choosing T* = <5 3a / 2 . 

In the isotropic case when R = R(\x\, k) we may simplify the above expressions for the diffusion 
matrices D mn and a mn . In that case we have 



^ mn(k, k) 



d 2 R{H' {k)sk,k) 

dXndXvn 



ds 



^HR"(H' (k)s,k)+[ 



k n k m \ R'(H' Q (k)s,k) 



k 2 



H^h)s 



ds 



R'(s,k) 



dS ( 5; 



h h 



H' (k) 

so that the matrix [D mn (k, k)] has the form 

D(k,k) =D (k) (l-k<g>k), D (k) 

In that case the functions \j are given explicitly by 

\H>(k)\ 2 \k\ 2 k 



k 2 



H' (k) 



R'(s,k) 



ds. 



and 



Q"nm{k) 



(d-l)D (k) 

K(k)\*\k\ 2 



T d ^(d-l)D (k) J sd 



D (k) 



k n k m d^l{\^) 



R'(s,k) 



ds 



\H' (k)\ 3 \k\ 2 
d(d - l)D (k)' 



2.6 A formal derivation of the momentum diffusion 

We now recall how the diffusion operator in 1)2.10(1 can be derived in a quick formal way. We represent 
the solution of l|2.8|) as 4> s (t, x, k) = if} S (5t, <5x, k) and write an asymptotic multiple scale expansion 
for ip s 

rP s (t, x, k) = fat, x, k) + y/Sfa (t, x, ~ k) + 5 fa (t, x, ~, k) + . . . (2.18) 
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We assume formally that the leading order term <j) is deterministic and independent of the fast 
variable z = x/<5. We insert this expansion into Q2.8[) and obtain in the order O (<5 -1 / 2 ): 



V z #i(z, k) • V k - H' (k)k ■ V z <f>i = 0. 



(2.19) 



Let 9 <C 1 be a small positive regularization parameter that will be later sent to zero, and consider 
a regularized version of (|2.19j) : 

1 



— — V z #i(z, k) • V k <j) - k • V z 0i + 00i = 0, 
H (k) 



Its solution is 



L (z,k) 



1 



#o(*0 



[ 0w m=1 
The next order equation becomes upon averaging 
d4> ^(dHxfak) 



yy 8Hi (z + sk, k) d<j)(t,-X.,\s) e ~0s di 



dz. 



8kr, 



(2.20) 



dt 



E 



V dk 



k • V z 0i - E (V z Fi(z, fc) • V k 0i) + #oW k • Vx^. (2.21) 



The first two terms on the right hand side above may be computed explicitly using expression Q2.2UJI 
for (pi : 

( 9Hl ^ k \ • Vzfa) - E (V z #i(z, fc) • V k 0i) 



dk 



-E 



A 8Hi{z,k) ~ 9 f 1 
/ ^ "to: 



5z m \H^(k) J 



dk 



dHi(z + sk, fc) 5(/)(t, x, k) -e s ^ 



+E 



E 



dHt(z,k) d ( 1 



9z n dk n 
d#i(z + sk, k) d<j>(t,x, k) 



Using spatial stationarity of H\(z, k) we may rewrite the above as 



e" 9s ds 



-E 



-E 



E 



9fTi(z,fc)f 9 / 1 

-/v« 



dHi(z + sk, fc) d(j)(t,x., k; _ es 



a a ( i 



dz 



dk n 



e^ s ds 



m,n=l 

E 

m,n=l 



<9z m <% m I #o(fc) 



dHi(z + sk, fc) <9^(i, x, k) _ 6)fj 



e" ys ds 



9 
dk m 





d 2 H 1 {z + sk, fc) ^ 00(t, x, k) fl 



dz n dz v 



dkn 



00 d 2 R(sk,k) 00(i,x,k) fl 



dk m \HUk)J dx n dx m dk n 

n,n=l \ u x ' 

d f 1 f°° d 2 R(sk,k) 00(i,x,k) 



1 E 

2 \ Hr.(kj J _ oc dx n dx m dk n 
m,n=l \ u 



ds , as 0-> 0+. 
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We insert the above expression into (|2.21j) and obtain 

% = E £(jWM)£)+fll(*)t.V,# P.22) 

m,n=l v 7 

with the diffusion matrix D(k, k) as in (|2.9j) . Observe that ()2.22|) is nothing but (|2.1()|) . However, 
the naive asymptotic expansion (|2.18|) may not be justified. The rigorous proof presented in the 
next section is based on a quite different method. 

3 From the Liouville equation to the momentum diffusion. Esti- 
mation of the convergence rates: proof of Theorem 12.11 

3.1 Outline of the proof 

The basic idea of the proof of Theorem 12.11 is a modification of that of p] |H] • We consider the 
trajectories corresponding to the Liouville equation ()2.8(l and introduce a stopping time, called t$, 
that, among others, prevents near self-intersection of trajectories. This fact ensures that until the 
stopping time occurs the particle is "exploring a new territory" and, thanks to the strong mixing 
properties of the medium, "memory effects" are lost. Therefore, roughly speaking, until the stopping 
time the process is approximately characterized by the Markov property. Furthermore, since the 
amplitude of the random Hamiltonian is not strong enough to destroy the continuity of its path, it 
becomes a diffusion in the limit, as 5 — > 0. We introduce also an augmented process that follows the 
trajectories of the Hamiltonian flow until the stopping time t$ and becomes a diffusion after t = t$. 
We show that the law of the augmented process is close to the law of a diffusion, see Proposition E31 
with an explicit error bound. We also prove that the stopping time tends to infinity as 5 — > 0, once 
again with the error bound that is proved in Theorem 13.61 The combination of these two results 
allows us to estimate the difference between the solutions of the Liouville and the diffusion equations 
in a rather straightforward manner (see Section 13.9)1 : they are close until the stopping time as the 
law of the diffusion is always close to that of the augmented process, while the latter coincides with 
the true process until r^. On the other hand, the fact that t$ — > oo as 5 — > shows that with a large 
probability the augmented process is close to the true process. This combination finishes the proof. 

3.2 The random characteristics corresponding to (12. 8fl 

Consider the motion of a particle governed by a Hamiltonian system of equations 

' «^V k F j )(^ im (*)(t;x,k)) 

< ^ = -^)(^,mM( f ;x,k)) (3.1) 
zW(0;x,k) = x, (0; x, k) = k, 

where the Hamiltonian Hs{x., k) := Ho(k) + ^/5Hi(~x., k), k = |k|. The trajectories of (|3.1|) are 
the characteristics of the Liouville equation (|2.8|) . The hypotheses made in Section imply that 
the trajectory (z( s \t; x, k), (t; x, k)) necessarily lies in C(T,5) for each T > 0, 5 G (0,5*(M)], 
provided that the initial data (x, k) E A(M). Indeed, it follows from the Hamiltonian structure of 
(|3.1|) that the Hamiltonian H$ (x, m) = Hq (to) + VdH± (z, to) must be conserved along the trajectory. 
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Hence, the definition J22J) implies that M^ 1 < |m^ (•; x, k)| < M s . We denote by Qf )X)k (-) the law 
over C of the process corresponding to (|3.1|) starting at t = s from (x, k) (this law is actually 
supported in C(5)). We shall omit writing the subscript s when it equals to 0. 



3.3 The stopping times 

We now define the stopping time rg, described in Section l3~Tl that prevents the trajectories of ()3.1j) 
to have near self-intersections (recall that the intent of the stopping time is to prevent any "memory 
effects" of the trajectories). As we have already mentioned, we will later show that the probability 
of the event [ rg < T] for a fixed T > goes to zero, as 5 — ► 0. 

Let < ei < £2 < 1/2, €3 G (0, 1/2 — €2), £4 £ (1/2, 1 — ei — €2) be small positive constants that 
will be further determined later and set 



N=[S~ ei ], p=[5~ e2 ], q = p[5- t % N 1 =Np[d~ 



(3.2) 



We will specify additional restrictions on the constants ej as the need for such constraints arises. 
However, the basic requirement is that ej, i = 1, 2, 3 should be sufficiently small and 64 is bigger than 
1/2, less than one and can be made as close to one as we would need it. It is important that e± < €2 



so that N <C p when 5 <C 1. We introduce the following {M. )(>o-stopping times. Let t k 
be a mesh of times, and ir £ C be a path. We define the "violent turn" stopping time 



(p) 



kp 



Ss(tt) : = inf 



t > : for some k > we have i £ 



,(p) ,(p) 



and 



(3.3) 



A^) • k{t) < 1 - 1 or A (t« - J-J • A(t) < 1 - I 
where by convention we set A(— 1/p) := A(0). Note that with the above choice of €4 we have 



A [t^ - 1/NA ■ K{t^>) >1-1/N, provided that 6 £ (0, 8 ] and S is sufficiently small. We adopt 

in ()3.3|) a customary convention that the infimum of an empty set equals +00. The stopping time 
S$ is triggered when the trajectory performs a sudden turn - this is undesirable as the trajectory 
may then return back to the region it has already visited and create correlations with the past. 
For each t > 0, we denote by Xt(iv) := [j X (s;ir) the trace of the spatial component of the 



0<s<t 



path 7r up to time t, and by Xt(q; vr) 
We introduce the stopping time 



[x : dist (x, Xt(ir)) < 1/q] a tubular region around the path. 



U 5 (tt) := inf 



t > : 3k >1 and t £ [*jf\ijf+l) for which G X t M ^ 



(3.4) 



It is associated with the return of the X component of the trajectory to the tube around its past - 
this is again an undesirable way to create correlations with the past. Finally, we set the stopping 
time 

r 5 (vr) := S s (tt) A U s (tt). (3.5) 



3.4 The cut-off functions and the corresponding dynamics 

Let M > be fixed and p, q, N, N\ be the positive integers defined in Section 13.31 We define now 
several auxiliary functions that will be used to introduce the cut-offs in the dynamics. These cut-offs 
will ensure that the particle moving under the modified dynamics will avoid self-intersections, will 
have no violent turns and the changes of its momentum will be under control. In addition, up to 
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the stopping time Tg the motion of the particle will coincide with the motion under the original 
Hamiltonian flow. 

Let at = 2 and a 2 = 3/2. The functions ipj : R d x Sf' 1 -> [0, 1], j = 1,2 are of C°° class and 
satisfy 



^■(k,l) 



' 1, if k ■ 1 > 1 - 1/JV and Mj 1 < |k| < M a 

(3.6) 

0, if k ■ 1 < 1 - aj/N, or |k| < (2M 5 )-\ or |k| > 2Mj. 



One can construct ipj in such a way that for arbitrary nonnegative integers to, n it is possible to find 
a constant Cmn for which HV'i'Hmn < C m . n iV m+n . The cut-off function 

®(t k-vr) ■= J ^ M ^ 2 M " ViVl )) fOT * G ^'ft) and ^ 1 (3 7) 

1 V 2 (k,^(0)) fort G [0,tS p) ) 

will allow us to control the direction of the particle motion over each interval of the partition as well 
as not to allow the trajectory to escape to the regions where the change of the size of the velocity 
can be uncontrollable. 

Let 4> : M. d x M. d — > [0, 1] be a function of the C°° class that satisfies (f>(y, x) = 1, when |y— x| > 3/q 
and (p(y, x) = 0, when |y — x| < 2/q. Again, in this case we can construct (j) in such a way 
that ||^>|| m ,n < Cq m+n for arbitrary integers m,n and a suitably chosen constant C. The function 
(f)k : W 1 x C — > [0, 1] for a fixed path it is given by 

Mr,*)= II ^( y ' x (^))- (3 ' 8) 

0<//g<4 P ii 

We set 

f 1, if < t < . . 

[ fc (y;vr), if *X_ < t < t^. 

The function $ shall be used to modify the dynamics of the particle in order to avoid a possibility 
of near self-intersections of its trajectory. 

For a given t > 0, (y, k) G R^ d and vr G C let us denote 0(t, y, k; vr) := k; vr)<I> (t, y; vr) . The 
following lemma can be verified by a direct calculation. 

Lemma 3.1 Let (/?i,/?2) be a multi-index with nonnegative integer valued components, m = + 
| /?2 1 - There exists a constant C depending only on to and M such that \d^d^e(t,y,k-7r)\ < 
CTl/3i| g 2|/3 1 | iV |/3 2 | for aUte [o ; T], (y,k) G „4(2M), vr G C. 

Finally, let us set 

F s (t,y,l;ir,u) = 0(t,5y,l;7r)V y fl'i(y,|l|;w). (3.10) 

For a fixed (x, k) G R^, 5 > and w G fi we consider the modified particle dynamics with the 
cut-off that is described by the stochastic process (y( s \t; x, k, u), (t; x, k, u>))t>o whose paths are 
the solutions of the following equation 



dt 



flS(|lW(t;x,k)|) + V5flkHi (^^,|ZW(i;x,k)|)l Z (dj (t;x,k 



= -±F S (t, «^^ |i W(t;x l k); V W(.;x,k) > lW(.;x,k)l (3-H) 



k y^(0;x,k) =x, 1^(0; x,k) = k. 
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We will denote by Q~\. the law of the modified process (y^ (•; x, k), x, k)) over C for a given 



5 > and by iS^k ^ ne corresponding expectation. We assume that the initial momentum k G A(M). 

(3.12) 



From the construction of the cut-offs we immediately conclude that 



3.5 Some consequences of the mixing assumption 

For any t > we denote by Tt the cr-algebra generated by (y^(s), l^ s \s)), s <t. Here we suppress, 
for the sake of abbreviation, writing the initial data in the notation of the trajectory. In this section 
we assume that M > is fixed, X\,X2 : (R x R rf x R^ 2 ) 2 — > R are certain continuous functions, Z 
is a random variable and fli, <?2 are R rf x [M _1 , M]-valued random vectors. We suppose further that 
Z,gi,g2, are .Fj-measurable, while X\,X2 are random fields of the form 



Xi(x,k) = Xi (^ J ff 1 (x,A : ),V x a^i(x,fc),V 2 c ^F 1 (x,fc) 



j=o,i 



For i = 1,2 we denote g { := (gf ] \gf>) where gf ] G R d and G [M _1 , M]. We also let 
C/(0i,0 2 ) :=E [Xi(0i)X 2 (0 2 )l , 0i,02 GR d x [M _1 ,M]. 



(2) 



(3.13) 



The following mixing lemma is useful in formalizing the "memory loss effect" and can be proved in 
the same way as Lemmas 5.2 and 5.3 of 



Lemma 3.2 (i) Assume that r, t > and 



inf 

u<t 



(i) y (5) {u) 
h s 



r 



(3.14) 



-a.s. on the set for i = 1,2. Then, we have 



E 



Xi( 9l )X. 



2U/2 



E[C/( 5 i )52 )Z] <2^(^j ||*i|U»raLoo||Z|| £ i ( n). 



(3.15) 



(ii) Let EXi(0, fc) = /or a// G [M , M]. Furthermore, we assume that 52 satisfies \3.1J^ , 



inf 



Si 



(i) y (s) (u) 



> 



r + n 



(3.16) 



and |<7^ — I > r i<5 1 / or some ri > 0, P-a.s. on i/ie event Then, we have 



E 



Xi( 5 i)X 2 ( 52 )z] -E[C/( 51 , 52 )^]| <C0 1/2 (^)^ 1/2 g) 11*2 IMI^IlLHn) (3-17) 



/or some absolute constant C > 0. #ere i/te function U is given by (|.V. J.'j]) . 
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3.6 The momentum diffusion 

Let k(i) be a diffusion, starting at k G at t = 0, with the generator of the form 



£F(k) = £ D m „(k,|k|)a fc 2 m , fcn F(k) + ^i? m (k,|k|)5 fem F(k) (3.18) 

m,n=l m=l 
d 

= E <^ (^,n(k, |k|)5 fen F(k)) , F G C °°(Mf). 

m,n=l 

Here the diffusion matrix is given by Q2.9|) and the drift vector is 

E rn(kl) = --J(KjJ2 I S d ^ S ^ ds, m = l,...,d. 



H'oW^Jo dx m dx 2 n 

Employing exactly the same argument as the one used in Section 4 of it can be easily seen that 
this diffusion is supported on §^ -1 , where k = |k|. Moreover, it is non-degenerate on the sphere, for 
instance, under the assumption ()2.6|) . cf. Proposition 4.3 of ibid. 

Let x ,k be the law of the process (x(t),k(£)) that starts at t = from (x, k) given by x(i) = 
x + J Q H' (\]^(s)\)^.(s)ds, where k(t) is the diffusion described by ()3.18j) . This process is a degenerate 
diffusion whose generator is given by 

£F(x,k) =£ k F(x,k)+^(|k|)k-V x F(x,k), F G C °°(Rf ). (3.19) 

Here the notation stresses that the operator C defined in ()3.18j) acts on the respective function 
in the k variable. We denote by 9Jl Xj k the expectation corresponding to the path measure £} x ,k- 

3.7 The augmented process 

The following construction of the augmentation of path measures has been carried out in Section 
6.1 of [Tl|- Let s > be fixed and n G C. Then, according to Lemma 6.1.1 of ibid, there exists 
a unique probability measure, that is denoted by 5 n ® s Qx(s),K(s), such that for any pair of events 
A G M s , B G M we have 5 n ® a Q X (s),K(s)[ A ] = M 71 ") and S v ® s Qx{s),K{s)[6s{B)] = 0. X ( s ),K(s)[ B ]- 
The following result is a direct consequence of Theorem 6.2.1 of |14j . 

Proposition 3.3 There exists a unique probability measure R^L on ^ suc -h that R^\\A\ := Q^^A] 
for all A G M T& and the regular conditional probability distribution of R^\ L [-\M TS ] is given by 
<5tt ®t 5 (tt) &x{T 6 (ir)),K(T S (n))> 7T G C. This measure shall be also denoted by Q^k ®n &x(t s ),k(t s )- 
Note that for any (x, k) G A(M) and A G M TS we have 

^liA] = Q%[A] = Qfl[A], (3.20) 

that is, the law of the augmented process coincides with that of the true process, and of the modified 
process with the cut-offs until the stopping time t$. Hence, according to the uniqueness part of 

Proposition E31 in such a case Q^L®^ &X(t 5 ),k(t s ) = Qxk®^ &x{t 6 ),k(t s )- We denote by the 
expectation with respect to the augmented measure described by the above proposition. Let also 
-^-ifkTT' E likn denote the respective conditional law and expectation obtained by conditioning i?^ 
on M TS . 
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The following proposition is of crucial importance for us, as it shows that the law of the augmented 
process is close to that of the momentum diffusion as 5 — > 0. To abbreviate the notation we let 

t 

N t {G) :=G(t,X(t),K(t)) -G(0,X(0),K(0)) - j (d e + C)G(q, X(g), K(g))) dg 



for any G G C 6 1,1,3 ([0, +oo) x R 2d ) and t > 0. 

Proposition 3.4 Suppose that (x, k) G A{M) and £ G Cb((M* d ) n ) is nonnegative. Let 70 G (0, 1/2) 
and let < ti < ■ ■ ■ < t n < < t < v < T . We assume further that v — t > J 70 . Then, there exist 
constants 71, C such that for any function G G C 1,l,3 ([T*,T] x M. 2d ) we have 

<L {Wv(G) - N t (G)] c}| < ct»(v - i)l|G|l5; 3 T] r 2 ^c (3.21) 

.Here C(7r) := £(X(ti), .RT(ii), . . . , X(i n ), K(t n )), tt G C(T, 5). T/ie choice of the constants 71, C does 
not depend on (x, k), 5 G (0, 1], C, fimes ii, . . . , t n ., T*,T, v, t, or the function G. 



Proof. Let = s < si < . . . < s n < t and B u . . . , B n G B(Rf) be Borel sets. We denote A := C 
and for any k G {1, . . . , n}, s < s k we define the events 

A fc := [tt : (X( S1 ),K( S1 )) G . . . , (X(* fc ), G 73 fc ] 

and their shifted counterparts 

4 S) := [tt : (X( Sfc - s ), K(s fc - s)) G B k , . . . , (X(s n - s),K{s n - s)) G B n ). 

For (x, k) G R 2d , vr G C and G G C 1 ' 1 ' 2 ^, +00) x M 2ci ) we let 

£ t G(t, x, k; tt) := H' (\k\) k ■ V x G(t, x, k) + Q 2 (t, X(t), K(t);n)jC k G(t, x, k) 

d 

-e(t,X{t),K(t);n) d Km e(t,X(t),K(ty,Tr)D m , n (k,\k\)d kn G(t,x,k) 

m,n=l 

and 

t 

N t (G) := G(t, X(t),K(t)) - G(0, X(0), K(0)) - J (d e + C e )G(g, X(g),K(g); tt) dg. 



It follows from the definition of the stopping time 75(71") and the cut-off function that 

V K Q(t,X(t),K(t);ir) =0, tG [0,t 5 (tt)], 

hence 

CtG(t,X(t),K(t)-,iT)=CG(t,X(t),K(ty,ir), t G [0, r^vr)]. 
We need the following result. 

Lemma 3.5 Suppose that (x, k) G »4(M) and £ G Gfe((]R 2d ) n ) is nonnegative. Let 7q G (0,1), 
0<ii<---<i n <T*<i<t;<T and t — T* > o~ 7 ° . Then, there exist constants j[, C > snc/i 
t/iat /or any function G G C 1,1,3 ([T*,T] x M 2d ) we have 

E%{[N V (G) - N t (G)]c}\ < &rt{v - t)\\Gff;J ] T 2 E^lc. (3.22) 

The choice of the constants 7J, C" does noi depend on (x, k), 5 G (0,1], iimes t\, . . . ,t n ,T^, T,v,t, 
or function G. 
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The proof of this lemma follows very closely the argument presented in Section 5.3 of and we 
postpone it until the Appendix. In the meantime we apply this result to conclude the proof of 
Proposition 13.41 We write 



n-l 



KlA N v(G) - N vArs(n) (G), A n ]=]Tl [Wp+l) (T 5 (7r))l Ap ^ 

p=0 

+ 1 K,i>)( r <s(7r))lA, l (7r)9Jt X ( 7 - 5 ( 7r )) iK(7 - (5 ( 7r ))[iV t ,_ 7 - (5 ( 7r) (G)]. 



>+i J 
(3.23) 



When t$(it) £ [s p , s p +i) we obviously have 

^X(t s {tt)),K(t s (tt))[ N v-t 6 {tt)(G), Aj^ ] = yXx(T S (TT)),K(T S (TT))[ N t-T S (n)(G) , A K p 

and ^x(t s (-k)),k(t s M)Wv-t s (-k)(G)] = 0. Hence the left hand side of ()3.23j) equals 

n-l 

^l{s p ,s p+1 )(n{n))l Ap (n)m X (T S (TT)),K(T S (n))[Nt- TS (n){G), )] ] 



(3.24) 



p=0 



= E^l n [N t (G) - N tATs(n) (G),A n ]. 
We conclude from (l3~23l . (EHil) that 

<U^( G )> A n ] = E^l n [N vArs{w) (G) + Nt(G) - N tArs{7T) (G),A n ] 



(3.25) 



^k,7r [^(«Arj (tt)) Vt (C) , A n ] 



and therefore 



^[JV W (G),^] = E^L [4?k,J^(«Ar 4W M(G), A 



(3.26) 



(S) 

x,k 



4!k,. [A r (,Ar 4W )V < (G),A n ] ,T 5 (tt) < t 



E *l,n [N (vATs(7r))vt (G),A n ] , t 5 (tt) > t 



The first term on the utmost right hand side of ()3.26|) equals E^.[Nt(G),A n , t$ <t], while the 
second one equals ^ k [Nr vATs \ vt (G) , B] . Here B := A n n [t§ > t] is an .Mf-measurable event. 
Suppose that 7q G (70 + 1/2, 1) and let L := [5 _7 o] be yet another mesh size parameter. We define 



and note that 



a := L-\([L(v A r 5 )] + 2) V ([Li] + 2)] 

[Lv]+2 



E^[NAG),B}= 

P=[Lt]+2 



N p/L (G), B,a=^ 



Li 



(3.27) 



Representing the event [a = p/L] as the difference of [a > p/L] and [a > (p + 1)/L] (note that 
[a > ([Lf] + 3)/L] = 0) and grouping the terms of the sum that correspond to the same index p we 
obtain that the right hand side of 1)3.27(1 equals 



E *l [N m+ 2)/L(G), B[ 



[Lv]+2 
p=[Lt]+2 



^ P +i/L(G ! )-iV p/i (G), B,<t> 



p+ 1 



(3.28) 
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Since the event B n [a > (p + 1)/L] is A^^ p 1 ^ //L -measurable, from Lemma 13.51 we conclude that 
the absolute value of each term appearing under the summation sign in (|3.28f) can be estimated by 

C'\\G\\ lx ^L- l Q [ X\B] which implies 



E% [N a (G), B] - E® [JV^+^-a (G), B] < CfS* \\G\$lpT* Q%[B] 



[Lv] + 1 - [Lt] 



L 



A direct calculation using formulas (|3.1|) allows us to conclude also that both \N a (G) — N, 



and \N, 



{[Lt\+2)L- 



E { H [N( V A T5 )vt(G), B] - E^ [N t (G), B] < E^ [N a (G) - N {vATs)yt (G), B 

[N a (G), B] - E® k [N m+S)L -i(G), B) + E^ [N m+2)L -i(G) - N t , B] 



(G) - N t (G)\ are estimated by C\\G\\ [ ^'P S^' 1 / 2 . Hence, (since 7^ > 1/2 + 



(vArg)Vt 



(G)\ 



7o J 



=.(*) 



+ 



< C5^ \\G\\f{fT 2 Q ( X\B] (v-t)y 6™ 



(3.29) 



for a certain constant C > and 71 := min[7Q — 70 — 1/2, 7^]. From (|3.26|) . (|3.29|) and the observation 
just below (|3,26|) . we obtain 

E%[N V (G) - N t (G),A n }\ < CS» M^T 2 R^^} („ - t) V 

for a certain constant C > and the conclusion of Proposition 13.41 follows. □ 

3.8 An estimate of the stopping time 

The purpose of this section is to prove the following estimate for -R^ k [t$ < T] . 

Theorem 3.6 Assume that the dimension d > 3. Then, one can choose 61,62,63,64 in such a way 
that there exist constants C, 7 > for which 



R xli T s < T \< CS^T, V5 e (0, 1], T > 1, (x,k) e A(M). 
Proof. We obviously have 

[t 8 <T] = [U 5 <ts,U s <T]U[S 5 < t s , S s < T] 



(3.30) 



(3.31) 



with the stopping times Sg and Us defined in ()3.3|) and (|3.4j) . Let us denote the first and second event 
appearing on the right hand side of ()3.31|) by A(6) and B(5) respectively. To show that (|3.31l) holds 
we prove that the -R^k probabilities of both events can be estimated by Cfr'T for some C, 7 > 0: 
see (EM) . (I3~4TH) and ' (ETUI . 

3.8.1 An estimate of R { ^ k [A(5)} 

The first step towards obtaining the desired estimate will be to replace the event A(5) whose definition 
involves a stopping time by an event C(5) whose definition depends only on deterministic times, see 
(|3.32|) below. Next we use the estimate Q3.21JI of Pr op osition 13 . 41 for an appropriately chosen function 
G to reduce the question of bounding the R^\. probability of A{5) by an easier problem of estimating 
its x ,k probability (£2 x ,k corresponds to a degenerate diffusion determined by (TJinH))- The latter 
is achieved by using bounds on heat kernels corresponding to hypoelliptic diffusions due to Kusuoka 
and Stroock. 
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We assume in this section to simplify the notation and without any loss of generality that 
h*(M) = 1. Note that then 



A(5) C A{5) :-- 



x[l)-x( l - 



< - : 1 < i < j < [Tq], \i - j\ > - 
q p 



(3.32) 



and thus 



x\l\-x( % - 



3 

< - 



: l<i<j<[Tq], |i-j|>^. (3.33) 



Suppose that : R d -> [0, 1] is a C 00 — regular function that satisfies /(x) = 1, if |x| < Ah/q 
and / (5) ( x ) = 0, if |x| > 5/q. We assume furthermore that i,j are positive integers such that 
(j - i)/q G [0, 1] and ||/ (<5) || 3 < 2g 3 . For any x G M d and i/q<t< j/q define 



G,(i,x,k;x ) :=9Jt x , k /^ (xK-t 



xo 



Obviously, we have 

d t Gj(t,x, k;x ) + £Gj(t,x,k;x ) = 0. 

Hence, using Proposition l3.4l with u = j/q and t = i/q (note that v — t> 1/p > 5 €2 and e2 G (0, 1/2)), 
we obtain that there exists 71 > such that 



(5) 



x,k 



/M(*(i)-«)-e J (i,*(i).jr(i);« 



(3.34) 



< c 



|G^ v ,sxo)||^' /9] r 2 ^, Vff 6(0,1]. 



According to ^3] Theorem 2.58, p. 53 we have 

||G^, v ;x )||^ /<z] < C||/W|| 3 < C</ 3 < CS- 3 &+*\ j G {0,...,[gT]}. (3.35) 

Hence combining (|3.34() and H3.35() we obtain that the left hand side of (|3.34() is less than, or equal 
to C £7i-3(e2+€3) f or a n ,5 e (q, 1]. Let now i = j - - so that 1 < i < z < j < [Tq]. We have 



(5) 
x,k 



X 



X 



< 



p 



< E. 



(S) 

x,k 



f (S) x 



X 



(3.36) 



E. 



(S) 
x,k 



E. 



(S) 
x,k 



/W (,(.) 



\y=x(i/q) 

According to (|3.34() and (|3.35f) we can estimate the utmost right hand side of (|3.36|) by 

sup j<m x , k .f (5) (x - yj : x,y G R d , k G A(2M)| + C 5^-^ €2+£ ^T 2 
To estimate the first term in (|3.37j) we use the following. 



(3.37) 
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Lemma 3.7 Let p,q be as in (|3.2j) . Then, there exist positive constants C\, C2 and C3 such that 
for all xj£ R d , k G 4(2M), ;£{!,..., [pT]}, (5 G (0, 1] we have 



Q 



X 



< 5 - 



C*2 

< Ci ( ^-j- + e" C3P ) . (3.38) 



We postpone the proof of the lemma for a moment in order to finish the estimate of ii[ k [i(i5)]. 
Using (|3,38|) we obtain that the expression in (|3.37|) can be estimated by 

d + e - C3P ^j + C7 57i-3(€2+ e3 ) r 2 < ^(d-C^+deg + eX p {-C 3 (T e2 } + C5 7i-3( e2 +e 3 ) r 2_ 

Hence, from ()3.33|) . we obtain that 

R xil A (S)} < [Tq? (Ci<5 (d - C2)£2+de3 + exp {-C 3 cT e2 } + C <5 71 - 3(e2+e3) T 2 ) (3.39) 

< CT 2 C < 5(rf-2-C 1 )e 2 +(d-2)63 + 5-2(^+^3) |_ C3( j- £2 j. +^71-5(62+63)^ < C(5 72 T 4 



for 72 := min[(d — 2 — Ci)e2 + (d — 2)53,71 — 5(e2 + £3)] > 0, provided that €2 + €3 < 71/5 and 
£2 G (0, (d — 2)63/(61 + 2 — d)). Here with no loss of generality we have assumed that C\ + 2 > d. 
Recall also that d > 3. Now suppose that 73 G (0, 72). Consider two cases: T 3 < <5 -73 and T 3 > <5~ 73 . 
In the first one, the utmost right hand side of (|3.39j) can be bound from above by C<5 72_73 X '. In the 
second we have a trivial bound of the left side by 5 73 / 3 T. We have proved therefore that 

R^l[A(5)} < C5?T (3.40) 

for some C, 7 > independent of 5 and T. 

The proof of Lemma 13.71 We prove this lemma by induction on j. First, we verify it for 
j = 1. Without any loss of generality we may suppose that k = . . . , k&) and kd > (4dMg) . 
Let Ann : -> K, m, n = 1, . . . , d - 1, E m : -> R, m = 1, . . . , d be given by 



D pq {\) := D pq (k^H,k~Wk 2 -l 2 ,k), E p (l) := E p (k~ l \,k~ l \f k 2 - l 2 ,k), 

when 1 G Z := [1 G B^ 1 : /c^VA; 2 - / 2 > (^Ma) -1 ], Z = |1|. These functions are C°° smooth and 
bounded together with all their derivatives. Note also that the matrix D = [l) mn ] is symmetric and 
D£ • £ > A |£| 2 for all £ G M^" 1 and a certain A > 0. The projection K(t) = (K^t), . . . , K d {t)) 
of the canonical path process (X(t; 7r), K(t; n)) considered over the probability space (C, M, H x ,k )) 
where kg := (1, Vk 2 — I 2 ), with 1 G Z, is a diffusion whose generator equals £, see (|3.18jl . It can be 
easily seen that (Ki(t), . . . , Kd-i(t)) t>0 , is then a diffusion starting at 1, whose generator M is of 
the form 

d-l d-l 

AAF(l,x) := J> 2 F(1) +^2a q (l)d lg F(l), F G Co°(R d-1 ), (3.41) 

p=l q=l 
where a q (l), q = 1, . . . , d — 1 are certain C°°-functions and 

d-l 

X p (l):=Y,D l p i 2 (l)d lq , P=l,...,d-1. 



q=l 



1/2 

The (d — 1) X (d — 1) matrix [D pq (1)] is non-degenerate when 1 £ Z. Let 



d-l 

MF(l, X ) := ]Tx 2 F(l,x) +X F(l,x), F G C °°(] 
P =i 



L d_1 x M d ) 
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where Xq is a C°°-smooth extension of the field 



k 



q=l 



9=1 



It can be shown, by the same type of argument as that given on pp. 122-123 of PP, that for each 
(x, 1), with 1 G Z, the linear space spanned at that point by the fields belonging to the Lie algebra 
generated by Xq, . . . ,Xd-i is of dimension 2d — 1. One can also guarantee that the extensions 
Xq, . . . , Xd-i satisfy the same condition. We shall denote the respective extension of N by the same 
symbol. 

Set lo := (fci, . . . , kd-i)- Let 1Zi , 72- x ,i be the path measures supported on C d ~ x and C d,d ~ l 
respectively that solve the martingale problems corresponding to the generators M and M with the 
respective initial conditions at t = given by lo and (x, lo). Let r(t,x — y, li,l2), t G (0, +oo), 
x, y G M. d , l\,h G IR^" 1 be the transition of probability density that corresponds to 7£ Xj i . Using 
Corollary 3.25 p. 22 of jl()j we have that for some constants C, m > 



r(t,y,k,l) <Cf 



Vy G R d , k, 1 6 t G (0,1]. 



(3.42) 



Denote by tz(tt) the exit time of a path it £ C d 1 from the set Z. For any 7r G C 6 ^ 1 we set also 
rz(vr) = t z (K(-,7t)). Let 5 : if" 1 -> Sf" 1 be given by 



5(1) :=(/i,..., Z d _i, Vk 2 -l 2 ), 1 = (Zi, • • • , G B*" 1 , Z := |1| 



and let 5 : C d ' d ~ x -> C be given by 5(vr)(t) := (X(i; vr), Soi^f; vr)), t > 0. For any A G M fz we have 
^x,l [£ _1 (^)] = X) 5(i o )[^]. Since the event [\X{l/p) - y| < 5/q] n [f z > is A4 fz -measurable 
we have 



x,k 



< 



< n 



X,l() 



5 1 

< -,T Z > ~ 

q p 



1 



(3.43) 



Here Co a denotes the volume of M d . To obtain the last inequality we have used (|3.42|) and an estimate 
for non-degenerate diffusions stating that 1Z\ [rz < 1/p] < Ce~ C3P for some constants C, C3 > 
depending only on d and Ao, see e.g. (2.1) p. 87 of |T1]. Inequality (|3.43|) implies easily (|3.38j) for 
j = 1 with C\ = m. To finish the induction argument assume that (|3.38|) holds for a certain j. We 
show that it holds for j + 1 with the same constants C\,C2 and C3 > 0. The latter follows easily 
from the Chapman-Kolmogorov equation, since 



x,k 



X 



3 + 1 



5' 
< - 

1. 



Q 



y.l 



X I J - 



5 

< - 



Q ( -,x,k,dy, dl 

.P 



induction assumpt. 

< Ci 



n c 2 



+ e 



-c 3P 



Q [ -,x,k,dy,d\ ) = C\ 



P — + e- C3P 
q a 



and the formula (|3.38j) for j + 1 follows. Here Q(t, x, k, •, •) is the transition of probability corre- 
sponding to the path measure x ,k- d 
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3.8.2 An estimate of R { ^[B(5)} 

We start with a simple observation concerning the Holder regularity of the K component of any 
path vr G B{5). Let us denote p := 2M 5 ~ 1 iV'- 1 / 2 and 



D :-- 



vr G C(T,<5) : |if(t) - K(s)\ > p for some k s.t. ^ < T and t G [t^.tj^J, s G [tj^,* 



f(p) +(p)i 



where Ms has been defined in (|2.7|1 and iV in (|3.2j) . Suppose that 7r G B(5), then we can find 
* G [4 P) '4+J' s e 4 P) ] for which ^(*) - K(s) < 1 - 1/iV. This, however, implies that 

l^)-^)| 2 >4l^(*)-^)| 2 >^' 
5 S 

thus tt £ D. Hence the desired estimate of R^\;[B(5)] follows from the following lemma. 

Lemma 3.8 Under the assumptions of Theorem 1.9. (A there exist C, 7 > such that 



Proof. We define the following events: 
Fx : 
F 2 : 
F 3 : 
F 4 : 



R^l[D] < CTS7, \/5 G (0, 1], T > 1, (x,k) G -4(M). 



lAVt) - K(s)\ > p for some s, t G [0, T], < t - s < -, t < t s 

P 
2 

\K{t) - K(s)\ > p for some s, t G [0, 21, < t - s < -, s > r 5 

V 

\K(t s ) - K(s)\ > J for some s G [0, 21, < r s - s < -, t 5 < T 
2 p 



\K(t s ) - K(t)\ > - for some t G [0, 31, < t - r 5 < - 
2 p 



(3.44) 



Observe that D C [j Fi. Note that Fi,F% are M T& -measurable, hence 



8=1 



R%[Fi\ = Q K ^[F i ], i = l,3. 
On the other hand for i = 2,4 we have 

R xi.i F i\ = J ^X{T 5 {%)),K(T S (-iT))[Fi,7r]Q%(dTv) 

where for a given ir G C 
Fa,-k '■■ 

Since all Fi, i = 1, 3 and -F^, i = 2,4, 7r G C are contained in the event 
F : 



\K(t) - K(s)\ > p for some s,t G [0, (T - r^vr)) A 0], < t- s < - 

p 

P 2 
#(0) - K(t)\ > - for some t G [0, (T - 73(71-)) A , < t < - 
2 p 



2 
- if(a)| > - for some s,t G 0, 2 1 , < t - s < - 
2 p 



(3.45) 
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(|3.44j) would follow if we show that there exist C > and 7 > for which 

Q%[F] < CW for all (x,k) G A(M) 



(3.46) 



and 



&*A F l < CT51 for a11 ( x ' k ) e ^(Mi). 



(3.47) 

The estimate (|3.47|) follows from elementary properties of diffusions, see e.g. (2.46) p. 47 of |13j . 
We carry on with the proof of Q3.46JI . The argument is analogous to the proof of Theorem 1.4.6 
of [II] . Let L be a multiple of p such that L := [5 _7 °], where 7q G (1/2, 1) is to be specified even 
further later on. Let also := k/L, k = 0, 1, . . .. We now define the stopping times t^tt) that 
determine the times at which the K component of the path ir performs k-th oscillation of size p/8. 
Let tq(tt) := and for any k > 

r fc+1 (vr) := inf [s^ > r fc (vr) : \K(s[ L) ) - K(r k (n))\ > 

with the convention that r n +i = +00 when r n = +00, or when the respective event is impossible. 
Let N# := min[n : r n +i > T] and 5* := min[r n — r n _i : n = 1, . . . ,N#]. Then, for a sufficiently 

small 5q and 5 G (0, 5o) we have F C [6* < 1/p] so we only need to estimate Q^fk probability of the 
latter event. 

Let / : R d -> [0,1] be a function of C£°(R d ) class such that /(0) = 1, when |k| < p/16 and 
/(k) = 0, when |k| > p/8. Let also /i(-) := /(• — 1) for any 1 G Note that according to Lemma l3.5l 
we can choose constants A p , C > 0, where C is independent of p, in such a way that A p < CT 2 p~ 3 
and the random sequence 



S' 1 



v 



E. 



(<5) 
x,k 



iV + 1 
L 



M N/L 



N 



(3.48) 



is a Q^ k -submartingale with respect to the filtration (jVl N / L } N>Q for all 1 with |1| G ((3M$) 1 , 3M$) 
provided that 5 is sufficiently small. We can decompose 



<5& 



< « 



5* < -, N # < [5- a ] 
P 



+ 



' x,k 



5* <-,N # > [5- 
P 



(S) 



n - n-i < - 

P. 



where a > is to be determined later. We will show that 



Q%[N # > HI < Ce 1 I 1 



T 



^1/2(61+62) ' 



(3.49) 



(3.50) 



and 



Q {8) 



r n+ i -r n < [L5 e 2 ]/L 



M 1 



< C5~<T 2 , 



(3.51) 



for < 7 < min[e 2 -3ei/2,7^-(l + ei)/2]. From (I3~1H1) . (I3~i31) (ETHUj) and (I33B we further conclude 
that 



{5) 

^x,k 



<s* < 



1 



< CT 2 5^~ a + Ce T [l 



(3.52) 
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for some C > 0, independent of 8 G (0, 1] and T > To, provided that we choose a E (l/2(ei + £2), 7). 
This is possible if min[e2 — 3ei/2,7 — (l+ei)/2] > (ei+€2)/2, which is true if we assume €2 > 10ei > 
and 1 > 7 > (l + e2)/2 + ei. Now, by the argument made after (j3.39j) we can always replace the first 
term on the right side of (|3.52|l by CT5 71 . We can also assume that the second term on the right hand 
side of 1)3. 52|) is less than or equal to CT5" 11 . This can be seen as follows. Let (3 := a — l/2(e\ + £2). 
The term in question is bounded by CexpjT— Ci<5 _/3 } with C\ := inf pg ( ,i] P 1 l°g (1 ~~ P/fy~ ■ 
For > 2T/C\ we get that expjT — C\5~P} is less than or equal to exp {— Ci5 _/3 /2}, while for 
5~P < 2T/C\ the left side of (|3.52fl is obviously less than 2T8@ /C\. In both cases we can find a 
bound as claimed. This proves Q3.46JI and hence the proof of Lemma [3.81 will be complete if we prove 
(l3~50|) and ffTST} . 

To this end, let Q^l. , 7r € C denote the family of the regular conditional probability distributions 



that corresponds to Q^fk I ' I -^ Tn ]- Then, there exists a M. Tn measurable, null Q^lc probability event 
Z such that for each it $l Z and each 1 £ if the random sequence 

Sjv,^=SU[o,iv/£]M7r)), iV>0 

is an (-^ iV ^ L ) A r>o submartingale under Q^Ictt- Let ^".^ := T n+i^(T n (Tr) + 2[L5 e ]/ L), where e G (0, 1) 
is a constant to be chosen later on. We can choose the event Z in such a way that 



Q^{ 7r [T n ^>r n (7r)] = l, Vvrg-Z. 



(3.53) 



Let S N . 



^N n ' then the submartingale property of ( Sjv,; 



N>0 



and (|3.53|) imply that 



(3.54) 



provided that 70 > (1 + ei)/2. The latter condition assures that p > C/{LyJ~8) so that K does not 
change by more than p during the time 1/L. In consequence of (|3.54|) we have 



E 



05) 

X,k,7T 



fK(r n (n)) K T„ j7r + 



+ 2A# > 1, 



as T nt7T — T n (ir) < 2[LS 6 ]/L. Since 

fK(r n (n)) \ K \T n ^ + 



fK(r n (n)) (K (T„ i7r )) 



< 



c 



Lp5 i n 



we obtain from (|3.55|) 



2^ > 41. [1 - /%wi Cr„ lir ))] - 



c 



so in particular 



C 



> E. 



(5) 

X,k,7T 



^x,k,- 



[L5 e 

I- fx (t„(tt)) (-^ ( r n+l)) ) T n+1 < T"n(7r) H p" 

[L5 f 



7"n+i < T n (jr) + 



We have shown, therefore, that 
[L5' 



Tn+l 



AT' 1 



CT 2 S e 
< ^ + 



C 



(3.55) 



(3.56) 



< C(J e - 3e i/ 2 r 2 +5 7 o- {1+ei)/2 ) < C8^T 2 (3.57) 
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for 71 < min[e — 3£i/2,7q — (1 + ei)/2] and some constant C > 0. We can always assume that 
T 2 cW 2 < 1. If otherwise, we can always write Qx — ^ 7 ^ 4 an d Q3.4fi(l follows. In particular, 
selecting e := (ei + £2)/2, one concludes from 1(3.57(1 that 



E^l[eM-(rn+i-r n )}\M Tn ] <e 

L5(ei+e 2 )/2 



-5(ei+e 2 )/2 



T n +1 - T n < 



L 



FT"! 



< 1 



T n +1 - Tn > 



\ L6 (e 1 +e 2 )/2 ] 



L 



M 1 



+ C(1- 

(j(ei+62)/2 



-5(n+^)/2 



(3.58) 



provided that 5 is sufficiently small. From (|3.58f) one concludes easily, see e.g. Lemma 1.4.5 p. 38 
of [HI, that (53Uj) holds. 

On the other hand, taking e = £2 in 1)3.57(1 we obtain 1)3.51(1 with < 7 < min[e2 — 3ei/2,7Q — 
(1 + ei)/2]. Hence the proof of Lemma 13.81 is now complete. □ 



3.9 The estimation of the convergence rate. The proof of Theorem 12.11 

Recall that 4>s, <t> satisfy 1(2.8(1 . 1(2. 1U[) . respectively, with the initial condition <J)q. We start with the 
following lemma. 

Lemma 3.9 Assume that <f>o satisfies the hypotheses formulated in Section WIX Then, 

d 

IMlJS < ||0o||o,o, E llMllK ^ H0o]|a,o- (3-59) 

i=l 

Furthermore, there exists a constant C > such that for all T > 1 

II^HK < C||0o||i,2- (3-60) 

In addition, for any nonnegative integer valued multi-index 7 = (a\, 02,03) satisfying I7I < 3 we 
have 

E ll^ iii;i2 , fci3 0|lK<CTH||0o||i,4, (3.61) 

*li*2,«3 = l 

Proof. The estimates (|3.59(l follow directly from differentiating (|2.10|) with respect to x. To obtain 
the estimates l(3.fi()(l and Ij3.fi 1(1 we note first that the application of the operator C to both sides of 
(j2.1(J() and the maximum principle leads to the estimate \\C(j)(t, x, -)\\ (A(M)) < I|£0o||l°°(.4(m)) f° r 
all t > 0, hence we conclude bound (j3.fi()() . 

In fact, thanks to already proven estimate (|3.59() we conclude that ||£^>(i,x, •)IIl°°(A(M)) < 
C 1 1 <^>o 1 1 1 , 2 for some constant C > and all (t, x) £ [0,+oo) x M d . Let Z be as in the proof of 
Lemmal3~71 Define 5 : Z x [M"\_M] -> A(M) as 5(1, A;) := (1, VA: 2 - Z 2 ), where Z = |1|. Let also 
VKl, fc) = 4> o 5(1, fc). We have (£ k 0) o 5(1, fc) = A^(l, fc), see ggp . The L p estimates for elliptic 
partial differential equations, see e.g. Theorem 9.13 p. 239 of [2j allow us to estimate 

\\Hw^{z) < c(U\\ LP(z) + W\\ LP{Z) ) < c||0ol|i,2. 

Choosing p sufficiently large we obtain that Yli ||^i^||z,°°(z) < C||0o||i,2> which in fact implies 
that ||D(-)Vk^(t, OlU^fSfz)) — C 1 1 00 1 1 1,2 - Obviously, one can find a covering of A(M) with charts 
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corresponding to different choices of the components of k being projected onto the hyperplane M. d 1 
and we obtain in that way that ||D(-)Vk0(i, — C 1 1 00 1 1 1,2 for all t > 0. Since the rank 

of the matrix D(k, k) equals d — 1, with the kernel spanned by the vector k, we obtain in that way 
the L°° estimates of directional derivatives in any direction perpendicular to k. We still need to 
obtain the L°° bound on the derivative in the direction k, denoted by d n := kid^ + • • • + kddk d - 
To that purpose we apply d n to both sides of (|2.1()j) and after a straightforward calculation we get 
dtd n cj) = Cd n (j) — 2C] i 4> + C\4> + H^(k)k ■ V x <fi, where 

d 



E 

m,n=l 



d_ 

dk n 



d k D mn (k, k) 



dk r 



Note that D(k, £;)k = implies that <9fcD(k, fc)k = hence ||£i0(i, OIU^f-AfM)) — C 1 1 00 1 1 1,2 - We 
already know that C^cfi and || Vx.^\\l oo (A(m)) are bounded, hence \\d n 4>(t, ■)||z, oo (.A(M)) < CH^olli^T" 
for t E [0, T]. We have shown therefore that \\4>(t, -)||i,i < C||0ol|i,2^ fo r t G [0>^1- The above 
procedure can be iterated in order to obtain the estimates of the suprema of derivatives of the 
higher order. □ 

Proof of Theorem 12.11 Let u £ [<5 7 o,T], where we assume that 7q (as in the statement of 
Lemma l3,5j) belongs to the interval (1/2, 1). Substituting for G(t, x, k) := (f)(u — t, x, k), £ = 1 into 
(|3.21[) we obtain (taking v = u, t = <5 7 o) 



E. 



{S) 
x,k 



<j>o(X(u),K(u)) - cf>(u- P°,X(P°),K(F°)) - / (d g + C Q )G{g,X{g),K{g))dg 



5 7 o 



< c W g \\i5 sj1t2 i v<5e(o,i]. 

Using the fact that |X(<5 7 o) - x| < C<5 7 o, \K(S^>) - k| < C<P° -1 / 2 , Qx k -a.s. for some deterministic 
constant C > 0, cf. (|3.11() . and Lemma 13.91 we obtain that there exist constants C, 7 > such that 



(3.62) 



E. 



(S) 
x,k 



MX(u), K(u)) - 4>{u, x, k) - / (d e + C e )G(g, X(g), K{g)) dg 



< C||G||^\j<FT 2 , 5 G (0,1], T > 1, u e [0,T]. 
We have however 



(3.63) 



4?k [0o(X(u),#(u)) - 0(«,x,k), r<5 > T] 



4 5 k [MX(u),K(u)) - 0(n,x,k), r 5 > T\ 



< C||G||S°$<TT a + (2||0 o ||o,o + T\\G\\^n [r a < T] 



[0,T]\ AW 



(3.64) 



Using /Vf 7 " 4 measurability of the event [t$ < T] we obtain that Q^k [ r <5 < ^1 = -^fL [ r 5 < ^1 an d by 
virtue of Theorem 13.61 we can estimate the right hand side of (ITMl) by 



Lemma 



C||G||^ 7 T 2 + C5 7 r(^2||^ || 0i o + T||G||^JJ < "' C(PT 5 . 

On the other hand, the expression under the absolute value on the utmost left hand side of (|3.64j) 
equals 



E. 



x,k 



b (X(u),K(u)) -0(u,x,k) 



E_ 



x,k 



J (X(«),Jf(u))-fcx ) k), r^T] 
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The second term can be estimated by 

m l3~30l 

\\<j>o\\o,oT, 

by virtue of Theorem 13.61 Since 

E^(|,|,k) =E0 o (zW(n;x,k),m^(n;x,k)) = E [ ^(X(u),K{u)) 

we conclude from the above that the left hand side of (|2.11|) can be estimated by C5 7 ||0o||i aT for 
some constants C, 7 > independent of 5 > 0, T > 1. The bound appearing on the right hand side 
of 1)2.11(1 can be now concluded by the same argument as the one used after (|3.39j) . □ 



4 Momentum diffusion to spatial diffusion: proof of Theorem 12.51 

We show in this section that solutions of the momentum diffusion equation (|2.1()|) in the long-time, 
large space limit converge to the solutions of the spatial diffusion equation (|2.12j) . We first recall the 
setup of Theorem 12.51 Let <^ 7 (t,x, k) = ^(t/j 2 , x/7, k), where </> satisfies ()2.10|) and let w(t, x, k) 
be the solution of the spatial diffusion equation (|2.12|) . In order to prove Theorem 12.51 we need to 
show that the re-scaled solution 7 (t, x, k) converges as 7 — > in the space C([0, T]; L°°(„4(M))) to 
w(t, x, k), so that 

\\w(t) - ^(t)\\ Lao{A{M)) <c( 7 T + 7 1 / 2 ) l^olko, 0<t<T. (4.1) 

Proof of Theorem 12. 51 The proof is quite standard. We present it for the sake of completeness 
and convenience to the reader. The function 7 is the unique C}' 1,2 ([0, +00), M 2d )-solution to 



2' 
7 ■ 



dt 

m,n=l 



7 (O,x, k) = o (x,k), 
see Remark 12.31 We represent <f)^ as 

= w + 7^1 + 7 2 ^2 + R- (4-3) 



Here w is the solution of the diffusion equation (|2.12|) . the correctors w\ and W2 will be constructed 
explicitly, and the remainder R will be shown to be small. The first corrector w\ is the unique 
solution of zero mean over each sphere of the equation 

E ( D mn(K k)^j = -H>(k)k • V X W. (4.4) 



m,n=l 



It has an explicit form 



Wl (t 1 x,k) = f; X ,-(k)^^fi (4.5) 

j=l 3 



with the functions Xj defined in (|2.14|) . The second order corrector W2 is the unique zero mean over 

m of the equation 



each sphere §^ 1 solution of the equation 

d 



m,n=l 
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Note that the expression on the right hand side of (|4.6|) is of zero mean since thanks to (j2.12j) and 
equality (|2.13|) we have 

din f 

H' {k)k-v xWl dn(k). 



dt T d _i j S d-i 

Equations (|4,4j) and (|4.6|) for various values of k = |k| are decoupled. As a consequence of this 
fact and the regularity properties for solutions of elliptic equations on a sphere we have that wi, W2 
belong to C([0, T]; L°°(A(M))). More explicitly, we may represent the function W2 as 



j,l=l J 



The functions ipj m (k) satisfy 



E ^ ( Dmn ^ k ^ ) = -^o(fc)%X/(k) +a iI (fc). (4.7) 



m,n=l 



A unique mean-zero, bounded solution of (|4.7|) exists according to the Fredholm alternative combined 
the regularity properties for solutions of 1)4. 7 Jl on each sphere With the above definitions of w, 

wi, W2, equation (|4.2[) for </> 7 implies that the remainder R in ()4.3|) satisfies 



ivv) rn 1 > / 



m,n=l 

We re-write this equation in the form 

d 

F'Cfc'it • V../? - — 

at 7 

i?(0,x, k) = 0o(x,k) - 0o(x,k) - (0,x, k) - 7^2(0, x, k), 

where / := — 7<9juii — 7 2 <9ju>2 — 7-^o(^)k ■ V x u>2- Here, as before, i? is understood as the unique 
solution to flUD that belongs to C fe lil,2 ([0, +00), M 2d ). We may split R = R± + R2 according to the 
initial data and forcing in the equation: Ri satisfies 

dRi l r w /f v,* „ „ 1 d / .* 



m,n=l 

2, 



dt 

m,n=l 

R(0, x, k) = -7^1(0, x, k) - 7^2(0, x, k) (4.9) 
and the initial time boundary layer corrector R2 satisfies 

f-i«*)£.VA-^ £ JL(iWt*)§£) =0 (4.10) 
i? 2 (0,x,k) =0 o (x,k)-^o(x,k). 

Using the probabilistic representation for the solution of (|4.9|) as well as the regularity of w\ and W2 
we obtain that 

\\Ri(t)\\L°°{A{M)) < C-yT, < t < T. (4.11) 
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To obtain the bound for R2 we consider ii^i, x, k) := i?2(7 3 / 2 i, x, k). This function satisfies 

' m,n=l x 7 

i$(0,x,k) =0 o (x,k)-0 o (x,k). 
We also define i? 7 ^ the solution of 



^ £ A^fc*)^ =0 (4.12) 



^(0,x,k) =0 o (x,k)-0 o (x,k). 

The uniform ellipticity of the right hand side of H4.12|) on each sphere implies, see e.g. Propo- 
sition 13.3, p. 55 of that the function satisfies the decay estimate on each sphere 

(7 7 (d-l)/4 (7 7 (<i-i)/4 
ll^2(*)lli°°(Srf-i) < f ( d -i)/2 ll^ollii^-i) < t(d _i)/2 ll^ollxoo(sj-i) ( 4 - 13 ) 



for t G [0, T] and, similarly, 



Cry (dr-l)/4, 

|V x i^(t)|| i0 o(s<-i) < t (d_i)/ 2 II ^0 II 1,0- 



Furthermore, the difference = i?2 ~~ -^2 satisfies 

^ - 7 1/2 ^(A:)k • V x g 7 - £ JL ^ mn(k , fc)|Q = 7 1/2 ^(^)k • V X B% (4.14) 

g 7 (0,x,k) =0. 

We conclude, using the probabilistic representation of the solution of Q4.14JI . that 

lk 7 (*)IU-(^)) < C7 1/2 t||0o||i,o 

and thus 

II^2(7 3/2 )IIl°°(.4(A-/)) < H-R2 (^lU^^M)) + lk 7 (l)IU°°(^(Af)) 

<c(7 (d - 1)/4 ||^o||o,o + 7 2 ||^o||i,o). 
The maximum principle for (|4.1()|) implies that we have the above estimate for all t > 7 s / 2 : 

\\R2(t)\\L°°(A(M)) < C (7 (d - 1)/4 ||^o||o,o + 7 1/2 ||V x ^o||i,o) , * > 7 3/2 - (4.15) 
Combining (|4.3|) . (|4. and (|4.15f) we conclude that 

\\w(t) - hWh-WM)) ^ C (7T + 7 ( ^ 1)/4 +7 1/2 ) ll<Mi,o, 7 3/2 < t < T, (4.16) 
and thus (|4.1|) follows, as d > 3. This finishes the proof of Theorem 12.51 □ 
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5 The spatial diffusion of wave energy 



In this section we consider an application of the previous results to the random geometrical optics 
regime of propagation of acoustic waves. We show that when the wave length is much shorter than 
the correlation length of the random medium, there exist temporal and spatial scales where the 
energy density of the wave undergoes the spatial diffusion. We start with the wave equation in 
dimension d > 3 

and assume that the wave speed has the form c(x) = cq + Vdci (x). Here cq > is the constant 
sound speed of the uniform background medium, while the small parameter 5 <C 1 measures the 
strength of the mean zero random perturbation c\. Rescaling the spatial and temporal variables 
x = x'/5 and t = t' jb we obtain (after dropping the primes) equation (|5. 1|) with a rapidly fluctuating 
wave speed 

Cj(x) = c + V5c! (|) . (5.2) 

It is convenient to rewrite ()5.1|) as the system of acoustic equations for the "pressure" p = (f>t/c and 
"acoustic velocity" u = — V</>: 

^ + V( C(5 (x)p) = (5.3) 
dv 

-^ + C(5 (x)V-u = 0. 

We will denote for brevity v = (u, p) £ M. d+1 and write (j5.3ft in the more general form of a first order 
linear symmetric hyperbolic system. To do so we introduce symmetric matrices As and D J defined 
by 

A s (x) =diag(l, 1, l,c«5(x)), and D J = ej (g> e d+1 + e d+1 ® ej, j = 1, . . . , d. (5.4) 

Here e. m G IR^" 1 " 1 is the standard orthonormal basis: (e m )^ = 5 m k- 

We consider the initial data for ()5.3|) as a mixture of states. Let S be a measure space equipped 
with a non-negative finite measure \x. A typical example is that the initial data is random, S is the 
state space and [i is the corresponding probability measure. We assume that for each parameter 
C £ S and e, 5 > the initial data is given by (0, x; £) := (— eV(^o(x), 1/cs(x)<Pq{x)) and vf(t, x; £) 
solves the system of equations 

^ +X>(x)I^ (^Wvf (x)) = 0. (5.5) 

3=1 

The set of initial data is assumed to form an e-oscillatory and compact at infinity family as e — > 0. 
By the above we mean that for any continuous, compactly supported function ip : R d -> R we have 

lim limsup / |99v|| 2 dk — > and lim limsup / [vf| 2 (ix — > 

R^+oo £ ^ 0+ J\k\>R/ e R->+oo £ ^ 0+ J| X |> R 

for a fixed realization £ £ 5 of the initial data and each 5 > 0. In the regime of geometric acoustics 
the scale e of oscillations of the initial data is much smaller than the correlation length 5 of the 
medium: e 5 <C 1. 
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The dispersion matrix for (|5.5j) is 

d d 
i$(x, k) = AsWkjDiAslx) = i ^ c s (x)kjD j = ipj(x) (k tg> e d+1 + e d+ i ® kj , (5.6) 
i=i i=i 

where k = X]j=i % e i- The self-adjoint matrix (— zPg ) has an eigenvalue Hq = of the multiplicity 
d — 1 , and two simple eigenvalues 

f4(x,k) = ±p 5 (x)|k|. (5.7) 

Its eigenvectors are 

b^= (k^,0) , m = l,...,d-l; b± = ^(j^±e d+l ^J , (5.8) 

where k^ 6 R d is the orthonormal basis of vectors orthogonal to k. 

The (d + 1) x (d + 1) Wigner matrix of a mixture of solutions of (|5.5j) is defined by 

W?(t > x,k) = ^L ? / / e^M^x-^jOvfftx+^Ocfy/iK)- (5-9) 

It is well-known, see [5l lllllT2*] . that for each fixed 5 > (and even without introduction of a mixture 
of states) when W*{t = 0) converges weakly in 5' (R d x R d ), as e -> 0, to 

Wo (x, k) = u° + (x, k)b + (k) ® b+ (k) + u° (x, k)b_ (k) ® b_ (k) . (5.10) 

then W*(t) converges weakly in S'(R d x R d ) to 

U 5 (t,x,k) =u^.(i,x,k)b+(k)®b + (k)+ui(t,x,k)b_(k)®b_(k). 

The scalar amplitudes u^± satisfy the Liouville equations: 



d t u s ± + v k ^4 • v x 4 - v x ^4 • v k 4 = o, 

«4(0,x,k) = t4(x,k). 



(5.11) 



These equations are of the form (|2.8|) . written in the macroscopic variables, with the Hamiltonian 
given by (|5.7|) . 

One may obtain an L 2 -error estimate for this convergence when a mixture of states is introduced, 
as in (|5.9|) . In order to make the scale separation e <C 5 <C 1 precise we define the set 

£„:={(£,*): |ln £ r 2 / 3+ ^<5<l}. 

The parameter /x is a fixed number in the interval (0,2/3). The following proposition has been 
proved in Theorem 3.2 of pQ, using straightforward if tedious asymptotic expansions. 

Proposition 5.1 Let the acoustic speed cs(x) be of the form and such that the Hamiltonian 

i/^(x) given by \5. 7| ) satisfies assumptions Qff.ffl . We assume that the Wigner transform satisfies 

W*(0, x, k) -> Wo(x, k) strong in L 2 (R d x M d ) as B (e, 5) -> 0. (5.12) 
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We also assume that the limit Wo G C 2 (M d x M. d ) with a support that satisfies 

supp Wo(x,k) C A(M) (5.13) 
for some M > 0. Moreover, we assume that the initial limit Wigner transform Wq is of the form 
W (x,k) = ^u°(x,k)n g (k), n ? (k) = b g (k)®b 9 (k). (5.14) 

q=± 

Let f/ <5 (t,x,k) = Up(t, x, k)IIp(k), where the functions Up satisfy the Liouville equations (|5. 
P =± 

Then there exists a constant C\ > that is independent of 5 so that 

||W*(i,x,k) - £/ 5 (t,x,k)|| 2 < C(5) (e||^ ||^e Clt/53/2 +e 2 ||^ ||H3e Clt/53/2 ) + \\w!(0) - W \\ 2 , 

(5.15) 

where C{5) is a rational function of 5 with deterministic coefficients that may depend on the constant 
M > in the bound 13$ on the support of Wq . 



The Liouville equations (|5.11j) are of the form (|2.8j) . Therefore, one may pass to the limit 8 — »• in 



(|5.11() using Theorem 12. II and conclude that Eii^_ converge to the respective solutions of 



m,n=l 

with the initial conditions as in 1)5.11(1 . Here the diffusion matrix D(k) = [D mn (k)] is given by 

D mn (k) = -- \° J ds, (5.17) 

where -R(x) is the correlation function of the random field ci(x): E [ci(z)ci(x + z)] = R(x). Fur- 
thermore, it follows from Theorem 12.71 that there exists ao > so that solutions of (|5.11|) with 
the initial data of the form n^_(0,x, k) = u±(5 a x,\s) with < a < ao, converge in the long time 
limit to the solutions of the spatial diffusion equation. More precisely, in that case the function 
u s (t, x,k) = u s + (t/5 2a , x/5 a , k) (and similarly for u s _) converges as 5 — > to w(t,x,k) - the solution 
of the spatial diffusion equation 

£ = f <wM^, (5.18) 
at ^ ax n dxm 

m,n=l 

w(0, x, k) = u° (x; k) := — !— / u° (x, k)dfi(k). 
with the diffusion matrix a mn given by: 

X m (k)dQ,(k), (5.19) 

1 d-l Js 1 *- 1 

and the functions \j above are the mean-zero solutions of 

m,n=l x ' 
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Theorems 12. 1| 12.51 and 12.71 allow us to make the passage to the limit e, 5, 7 — ► rigorous. The 
assumption that e <C 5 <C 7 is formalized as follows. We let 

/C M)P :={m, 7 ) : 5> |lne|- 2 / 3+ ^ and 7 > , 

with < p < 2/3, p E (0, 1). Suppose also that E C^(R^) and supp Uq C A(M). Let 

VF°(x,k) := u^(x,k)b+(k) ®b+(k) + n° (x,k)b_(k) <g> b_(k), (5.21) 

and 

W{t, x, k) := w+(t, x; k)b + (k) <g) b+(k) + w;_ (t, x; k)b_ (k) <g) b_ (k). (5.22) 
Our main result regarding the diffusion of wave energy can be stated as follows. 

Theorem 5.2 Assume that the dimension d > 3 and M > 1 are fixed. Suppose for some < p < 
2/3, p £ (0, 1) we /iai>e, u>ii/i as in (pTgTj) and W/ denned by (Q 



2 

dxdk — ► 0, as (e, 5, 7) — ► and (e, (5, 7) E /C MiP . 



Then, there exists p\ E (0, p] suc/i that for any T > T* > we have 



sup 

te[T„,T] 



EW E M^,^k)-W(t, X) k 12 



dxdk — ► 0, as (e, (5, 7) — > and (e, 5, 7) E /C At 



in ■ 



Here W(t, x, k) is 0/ i/ie /orm with the functions w± that satisfy (|5. 18|) £/te initial data 

w± (0, x, k) = u± (x, k) . 

The proof follows immediately from Theorems 12.11 l2~5l and l2~7l as well as Proposition 15.11 

A The proof of Lemma 13. 5L 

Given s > a > 0, ir £ C we define the linear approximation of the trajectory 

L(a, s; tt) := X(a) + (s - a)H' (K(a))K(a) (A.l) 

and for any v E [0, 1] let 

R{v, a, s; vr) := (1 - v)L(a, s; vr) + vX(s). (A.2) 
The following simple fact can be verified by a direct calculation, see Lemma 5.4 of pQ. 
Proposition A.l Suppose that s > a > and tt E C(5). Then, 

s 

\X(s)-L(a,s;n)\ < D(2M s )y/6(a - a) + J \H' Q {K(p))k(p) - H' Q {K{a))K{a)\dp. 



31 



We obtain from Proposition IA.1I for each s > a an error for the first-order approximation of the 
trajectory 

\ z W{s)-lW{a,s)\ < D{2M s )VS(s - a) + C ( s 6 € (0,6*(M)]. 

2y 5 

Here l& (a,s) :=z^(a) + (s — a)rhS '(a) is the linear approximation between the times a and s and 

C:= sup (M s h* (M s ) + h* (M s ))D(2Ms). 

<5e(0,<5.(M)] 

With no loss of generality we may assume that x = and that there exists k such that t,u € 
[t^ , £ ) . We shall omit the initial condition in the notation of the solution to 1)3.11(1 . Throughout 
this argument we use Proposition IA. II with 

o-(s) :=s-^ 1_7A for some 7aG (0,1/16 A (l-e 4 )), se[t,u]. (A.3) 

The aforementioned proposition proves that for this choice of a we have 

\L^ 5 \a,s) -y {5) {s)\ < C A 6 3/2 ~ 2 ^ A , V5e(0,l]. (A.4) 

Throughout this section we denote C = ((y^ (h), . . . , y^ (t n ) , Z (<5) (*„)) . We assume first 

that G G C 2 (R^) and ||G]| 2 < +oo. Note that 



d u / 

G(l^(u)) - G(Z«(t)) = -^1LJ d 3 G{l {S \s))F h5 \s, 
We can rewrite then (|A.5|) in the form + + 1^ , where 
:=-^E/^)(,))F, 



yWis 



&-,l®{s) | ds. (A.5) 



i.,^,i<V)U 



/ {3) := ^ £ // %G(jW(p))^ ( fl ,»^« J«(p)) ^ (p,^^>(p)) e^p 

and cr is given by l)A.3|) . Each of these terms will be estimated separately below. 
A.l The term E[I^Q 

The term 1^ can be rewritten in the form + , where 

:= -^£/^(Z(V))F, 5 ^^p^,l^(a^j ds, 



and 

u 1 

j(2) . 

; J =1 t b 



*>J=1 t n V / 



(A.6) 
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where, see (EZQ) and D, L( 5 )(<7, s) = L(<7, s; y^(-) , I® (•)) , s) = R(<7, s; (■) , (■)) . We 

use part (i) of Lemma l3~2l to handle the term E[JW£]. Let -Xi(x, k) = —d Xi H\{x., k), ^(x, k) = 1, 



Z = Q\t 



f ) ,^ I i w w)a,G(iW(a))c 



and <7i = (L^ s \a, s)5 1 , |Z^(o")|). Note that 51 and Z are both measurable. We need to verify 
(SHU). Suppose therefore that Z + 0. For p G [0,i£ x ] we have |Z (<5) (cr, s) - > (2g)~\ 

provided that Ca5 3 ^ 2 ~ 2 " /a < l/(2q), which holds for sufficiently small 5, since our assumptions on 
the exponents €2,63,7,4 (namely that £2,63 < 1/8, ja < 1/8) guarantee that €2 + £3 < 3/4 — ta/2. 
For p G a ] we riave 



(£<V, s) - yW(p)) • (tfcU > (, - er)flS(|I«(a)|) Z W (a) • 2 W (ft 



(S) 



W (,{p) 



(A.7) 



+ 



2 



(5) (pi; 



> (s -a)K(2M 5 ) ( 1 



iV 



+ 



K(2M S ) - VSD(2M 5 )] (s-p)(l 



N 



> (s - a)K(2M 5 ) 1 



N 



provided that 5 G (0, So] and So is sufficiently small. We see from (jA.7j) that (|3.14[) is satisfied with 
r = (1 - 2/N) h Jf (2M s )S 1 ~^ A . Using Lemma l3~2l we estimate 



E[J«C] 



< 



g(2Mg) 



||G||iE[C] / </> C 



s — CJ 



,(1) 



r/.s 



(A.8) 



< cfwhm^'^^h-^) (u-t)< cfwGhnmu-t) 

and exists by virtue of assumption (|2.4|) . On the other hand, the term j( 2 ) defined by (|A.6 
be written as J^ 2 ) = j[ 2 ^ + , where 



may 



J 



(2) 



i £ 1 d 3 G^\a))d yi F h& U^p^^\a)\ (y?\s) - L?\a, s)) ds 

i,j=l t V / 



and 



J 



(2) 



1 ^ 
'55/2 



Mil 



(A.9) 







x djG(l® (a))(yf } ( S ) - Lf } (a, S ))(yf (s) - hf (a, s)) ds dv M. 



The term involving may be handled easily with the help of QA.4JI and Lemma 13.11 We have 
then 

|E[jf } C]| < C A 4) D(2M s )E[C}\\G\\ l (u - t )S^ /2 S^ A - 4{ - e2+e3) T 2 (AJ0) 
< C A 5) 5 1/2 ~^ A - i(e2+e3) T 2 (u - t)E[C]\\G\\x. 
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(2) (2) (2) (2) 

In order to estimate the term corresponding to J{ we write J\ = J\{ + J\ 2 > where 



7(2) .. 

J l,l ■- 



d u s / (g\ 



(T 



(A.ll) 



x( a -pi)^-[flS(l' (a) (pi)l)^(pi) 



ds dpi 



and 



42 — JdMWWM* (s,^^^Ha^J 8 lHl (^,\l(%^ lf\p)dsdp, 



with 



A 

dpi V 

+H' (\l^(p 1 )\)\l^( Pl )\- 1 



H^(p 1 )\)l\ S \p 1 )} =H'^(pi)\)(f\pi),^-l i5) (pi))^i ( i S) (pi) (A.12) 



d 



dpi 

^i?\pi)-f\pi)±i( s H P i)) R 4 5 \pi) 



d 
dpi 



We deal with j[ 2 ^ first. It may be split as J^l = J121 + ^122 + ^12 3' where 



d u 
m =1 * 



8 \ s, 



(A.13) 



It S 1 



(2) ._ 
1,2,2 •" 



'J^ 1 t o- ^ 

'R®(v,v,p) 



5 2 



{yf\p)-Lf\a,p))if\p)dsdpdv 



and 



(2) .. 
1,2,3 ■" 



A U S P 



i,i=l 




5 



t <J <J 



dpi 



ds dp dpi . 



By virtue of (|A.4[) . definition 1)3.10}) and Lemma 13, II we obtain easily that 

|E[jg )2 C]| < cf5 1/2 -^ A+2€2+2e3) \\G\\iT(u-t)EC 



(A.14) 



The same argument and equality QA.12|) also allow us to estimate |E[j{ 2 2 3 C]| by the right hand side 
of (EH- 
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Using Lemma 13.11 and the definition ()3.1U|) we conclude that there exists a constant C A > 
independent of 5 such that 

< C^d 1 - 2 ^^, i, j = 1, . . . , d. 
Therefore, we can write 

, d u 3 



(A.15) 



t a 



4wJri (^,|,«V)|W^,^(«r) 



)l if'MC 



(is dp 



< ci 8) 5 1-7A-2(6a+e3) (u - t)||G||iTEC 



With our choice of the exponents we have 5 < (2q)^ 7A 2 ( e 2+f3)] 1 f or a \\ $ g £ ) where <5 > is 
sufficiently small. We apply now part (ii) of Lemma 13.21 with 



z = a,G(i«(a))e (#U«(a,'),i (a V)) ifV)c, 

Xi(x, fc) := d 2 H x {x } k), l 2 (x) := 3 fc #i(x, fc), 



5i - 



52 := 



We conclude that 



where 



E 



r = 4 9) (p-a), n =4 9) ( S -p). 



J SV] + ^ E / / E L W ( ») (A.16) 



t <7 



Xd Vi,Vj R l 



LW(a, S )-£( 5 V,p) 



, |I ( V)I 1 f (*)< 



ds dp 



< Cf 5 l -"< A - 2 ^ +ei) {u - t)||G||iTEC 



r (io) } } 

+^-||G|| 1 E[C] J J ^ 2 

t a 



25 



6 1 ' 2 



Cf{s-P) 
25 



ds dp, 



i*i(y, k) := E[#i(y, k^H^O, k)], (y, k) G R d x [0, +oo). (A.17) 

We can use assumption 1)2. 4[) to estimate the second term on the right hand side of ()A.16|) e.g. by 
C { 2 1] 5(u - t)||G||iEC. The second term appearing on the left hand side of (|A.16|) equals to 



W E {djG(l W (<r))9 (4 P) , (a, I« (a) 

3=1 t 



(A.18) 
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s-p 



H> (\lW(*)\)i {S) (a),\lW(a)\) dp 



C> ds 



and integrating over dp we obtain that it equals 



(A.19) 



+E / E ( f^^ e (4 p U ( V,*m ( V)) ' ( V)> c} 

t 



By virtue of (j2.5j) the second term appearing in l)A.19j) is bounded e.g. by 5{u — t)\\G\\i¥,C, for 

( 12) 

some constant C \ > 0, thus we have shown that 



J'=l 



^(I* (< V)I) 



(A.20) 



< C , ^ 13) <5 1 -^- 2 ( e2+£3 )(« - t)||G||iTEC- 



Let us consider the term corresponding to cf. (|A.11|> . Note that according to (|A.12j) and 
1 1|) we have \ = jf'l 1 + 2 , where 



o" J 



with 



r, (p, y, 1) := |irX(|l|) 1, F 5 (p, y, 1) L - ify (p, y, 1) - i^(|l|) 1, F 5 (p, y, 1) _ I 



while 



J 



(2) ._ 
1,1,2 - § 2 



hj— 1 t a a ^ 



x {s-p^Yi (p u V^p±jV)(p 2 )\ dsd Pl dp 2 . 



(A.21) 



Note that l^r^l < C^ 14 ^5 1//2 for some constant > 0. A straightforward computation, using 

IpPJ) and Lemma mi shows that lE^g^C] | < C^ 15) J 1 / 2 -( 3 ^ +2e2+2e3 )(u - t)||G||iTE[C]. An appli- 
cation of (|A.4|) , in the same fashion as it was done in the calculations concerning the terms 2 C] 
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and E[J-p2 3 £], yields that 



d 



ij =1 t a 



I pi> — ^ — '* y a ) 







H 


ds dpi 



djGil^ia^dy^s (s LS( f ,S \ l^(a)) (A.22) 



< c { ^ ) 5 1,2 ^ A+2 ^ +2 ^\u - t)||G||iTE[C]. 



For j = 1, . . . , d we let 



i,k=l 
d 



+ E H 'o(m\- 1 dl yi , yi R(y-y',\l\), 

i=l 

and also 

A(t,y,y',U) :=G(^y,l;^)G(t,y / ,l;vr), t > 0, y, y' £ R d , 1 £ R* tt £ C, (A.23) 
P: = (l^^J^PiJ^V))^- (j-^^^-^^PiM^V)) and 

9» :=e( S ,^( S ),zW( S ); 2 /( 5 )(-),i {5) (-))- 

Applying Lemma 13.11 and part ii) of Lemma 13, 2\ as in (jA.15|) and (|A,16|) , we conclude that the 
difference between the second term on the left hand side of ()A.22|) and 



5 2 



E / / ( s ~ Pi) E \d 3 G(l^(a))A(a,P)V 3 (P s ) ( 



ds dp 



J =1 i <7 



can be estimated by C^ 1 ' ' 5"<a' (u — i)||G||iE[£] for some 7^ > 0. Using the fact that 
estimate (|A.4|) and Lemma 13. II we can argue that 



(A.24) 



(A.25) 



A(<7,P) - 8 (s) < ^(J 1 / 2 - 7 - 4 -^ + jV2-2(7A+62+e 3 ) r )_ 

We conclude therefore that the magnitude of the difference between the expression in (lA~24l) and 



d 



^G(Z(' 5 )(a))G i ( S ) /( S -^-(P,)dpi c 



(is, 



(A.26) 



can be estimated by C^ 19 ^ 7 ^ {u— t)\\G\\\T¥.[Q for some 7^ > 0. Using shorthand notation Q(a) :- 



H' (\l^(a)\)f\a) 



1 



we can write the integral from a to s appearing above as being equal to 

d 

E teai ( V)i) - p^(i^V)i)) O 



s-S 1 —*-* 



i,fc=l 



( s "Pi) 

xlf (a)lf (a)+H^\a)\)\l^(a)\- 1 J2dy i , ym R 



s - pi 



Q(aU {S \a)\ 



i=l 



S ~ Pi 



Q(a),\l {S \a)\ 



dpi, 
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which upon the change of variables p\ := [s — pi)/ 8 is equal to 

8—>a 



Pi 



o 



,fe=i 



i=l 



(A.27) 



Using the fact that 
d 



k=l 



H' (\l^(a)\) d Pl 



d 2 my .R( Pl Q{<T)M s Xa)\ 



we obtain, upon integrating by parts in the first term on the right hand side of (|A.27p . that this 
expression equals 



i=i 



H' Q {\l^{a)\r l [K{\l^{a)\) - H' (\l^(a)\)) ]T 
dl >y .R(p 1 Q(a),\l^(a)\)if\a)dp 1 



S-^d 2 yi>y .R(6-^Q(a),\l^(a)\)il 5 \a) 



<5"TA 



(A.28) 



<5"TA 



+^(|zW(a)|)|zW(a)r 1 / Pl dl y ^R( Pl Q(a),\l^(a)\) d Pl . 



(A.29) 



Note that V J R(0) = and 

Y,9l yj R( P iQ(a),\l^(a)\)i\ 5 \a) 



d 



H>(\l^(a)\) d Pl ; 
We obtain therefore that the expression in ()A.28|) equals 

d 



—d y .R( Pl Q(a),\l^(a)\ 



H'oQl^ia)])- 1 (HZ(\l( s \a)\)-H' (\l®(a)\) 



^5-^dl yj R («T^ Q(a), |Z( 5 )(a)|) zfV 



i=l 



-H' (\l( s Xa)\)-%.R(8-^Q(*),\l( s \a 



(A.30) 



i=l 



^(I^DI^V^E / Pidl, yi , Vj R ( Pl Q(a),\l^(a)\ ) dpi 



Recalling assumption (|2.5j) we conclude that the expressions corresponding to the first two terms 



(3) 

appearing in (|A.30|) are of order of magnitude 0(<5 7 a ) for some 7^ > 0. Summarizing work done 



„(3) 



38 



in this section, we have shown that 



E 



< cJfV^u - t)||G||iT^EC (A.31) 



< 4 > 



for some constants C^f°\7^ > and (cf. (|A.17j) ) 

CjQ) :=E J (l,\l\) + 



dy 3 Rl (0,|1|) 



d +oc 

Ej (l,k):=-^P-J2 I p 1 ^ y .R{p l H' Q {k%k) dp x , j = l,...,d. 

8=1 n 



A. 1.1 The terms Ef/^C] and E[I^(] 

The calculations concerning these terms essentially follow the respective steps performed in the 
previous section so we only highlight their main points. First, we note that the difference between 
E[/( 2 )<C] and 



1 d 



E 



t a 



dsdp (A.32) 



is less than, or equal to cf l] *Pa' ( u - t)\\G\\iE[(], cf. ifPSJl . Next we note that (TQ2l) equals 

u s 



„( 5 ) 



(5) 



w c 



ds dp 



1 d 
+ T2 E 



« s 1 
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E 



ij,fe=l t a 



s 

u s 1 



5 



(A.33) 



+ 



Y2 E /// E ^)(.))^F,J S ,^M,^(,) 



*J'.*=1 t (T 



x dy k Fi,s \ p, - (yf (p)-4 5) (a,p))C 



(is dpdv. 



A straightforward argument using Lemma l3.1l and ()A.4|) shows that both the second and third terms 

of (TQ3ll can be estimated by C , J 3) <y 1 / 2 -( 3 T- A+2<a+26s )(ti-t)||G||ir 2 E[C]. The first term, on the other 
hand, can be handled with the help of part ii) of Lemma 13.21 in the same fashion as we have dealt 
with the term jj 2 ^ 1 ; given by (lA~13l) of Section and we obtain that 



E 



d " 
3=1 * 



(A.34) 
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<4 24) ^ 6) (n-t)||G||iTE[C]. 



Here 



Dj(l):= Wo,D 



R 2 (y,l) -EldiH^y, 1)^(0,1)}, 

d 



(A.35) 



i=i 



e,( S ) := ^e( S ,yW( S ),iW( S );yW(0,i (<5) (-)). 
Finally, concerning the limit of E[/( 3 )£], another application of (|A.4|) yields 



E[I^(] —X < ci 25) (T^"(«-t)||G||iE[C] 



(T) 



where 



x:=i / / E {a^) ( .))F, 5 (p,^^,^(.))c} 



(A.36) 



Then, we can use part ii) of Lemma 13.21 in order to obtain 



d u 

J " E / A,i(i (5) (^),|i w WI)0 2 ( S )^.G(zW(a))d s 



< cfV^fa _ t )||G|| 2 TE[C]. (A.37) 



Next we replace the argument <r, in formulas (EOT!) . (lA~34l) and (|A~36l) . by s. This can be done 
thanks to estimate (jA.25|) and the assumption on the regularity of the random field H\(-, •). In order 
to make this approximation work we will be forced to use the third derivative of G(-). 
Finally (cf. (TOTl) . ifPSJl ) note that 

V y R! (0, + V y R 2 (0, 1) = V y | y=0 E® #i (y, (y , 0] = 

Hence we conclude that the assertion of the lemma holds for any function G G C 3 (Rf) satisfying 
1 1 <j 1 1 3 < +oo. Generalization to an ar bitrary G G C b M ' 3 ([0, +oo) x R^) is fairly standard. Let r be 
any integer and consider s k := t + kr^ l {u — t), k = 0, . . . , rp. Then 



E 



G(u,y^(u),l {5) (u)) - G(i,yW(t),zW(0)]C} 
^2E{[G(s k+1 ,y^( Sk+1 ),l^( Sk+1 )) - G(s k ,y^(s k ),l^(s k ))]c} 

ro — 1 

^E{[G( Sfc ,yW(^,Z (5) (^))-G( Sfc ,yW(^0,Z (5) (^))]C} 



(A.38) 



ro— 1 



fc=0 
ro — 1 



fc=0 
ro-1 



+ Y J ^{[G(s k+ i,y {5 Hs k+1 ),l^Hs k )) - G( Sk ,y^\ Sk )J^(s k ))]c} 



k=0 

Using the already proven part of the lemma we obtain 

r -l 



^E{[j? 5H1 (G( Sfcl »W( Sfc ) ) .))-^(G( S|fcl »W(« fc ) ) 0)]c} 



fc=0 



(A.39) 
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<C% r) 8^\u-t)\\G\\ 1X3 T 2 EC 
On the other hand, the second term on the right hand side of ()A.38|) equals 



rp— l 
k=0 



r»-V y (A.40) 



G(p,y®(p),l( s \s k ))(dp) 



The conclusion of the lemma for an arbitrary function G £ C. ' ' ([0, +00) x M% d ) is an easy conse- 
quence of (|A.38|) - (|A.40|) upon passing to the limit with ro — > +00. □ 
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